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SECTION  I 


nrCKODUCTIOM 


Philosophically,  aircraft  dynamic  flying  qualities  requirements 
should  be  Independent  of  the  method  of  control  system  mechanization* 
However,  from  a  practical  viewpoint,  It  Is  known  that  the  type  of 
control  mechanization  enq>loyed  does  affect  the  aircraft  response  char¬ 
acteristics.  Consequently,  there  Is  a  need  to  practically  assess,  via 
in-flight  evaluation,  the  extent  of  the  control  response  modifications 
Introduced  by  digital  Implementation. 

All  of  the  significant  digital  controller  implementations  (e.g., 
F-18  and  Shuttle)  as  well  as  the  system  presently  being  developed  for 
the  F-5G,  have  utilized  the  emulation  approach.  Emulation  Is  a  process 
wherein  an  existing  analog  control  law  Is  "digitized"  by  the  use  of  a 
so-called  "substitution-for-s"  rule,  usually  the  Tustln  transform.  The 
verisimilitude  of  this  approach  Is  critically  dependent  on  open-loop 
system  bandwidth  and  sample  rate. 

In  contrast,  it  is  the  objective  of  the  present  study  to  use  direct 
digital  design  procedures  and  synthesize  low  data  rate  control  laws  for 
a  conventional  coupler  (zero-order  hold)  and  a  more  complex  coupler 
which  introduces  less  control  roughness  (slewer).  These  control  laws 
will  then  be  flight  tested  to  evaluate  the  modifications  Introduced  In  a 
particular  direct  force  mode. 

Specifically,  the  analog  fly-by-wire  system  of  the  Variable  Sta¬ 
bility  Navlon  will  be  programmed  to  simulate  an  open-loop,  lateral- 
directional  YF-16  at  M  -  0; 8  and  an  altitude  of  20,000  ft.  A 
microprocessor-implemented  controller  is  then  closed  around  the  bare- 
alrframe  YF-16  simulation  to  force  a  wings-level  turn  (WLT)  mode.  These 
flight  evaluations  will  replicate  a  previous,  completely  analog, 
in-flight  simulation  which  will  serve  as  a  baseline  configuration. 


An  ultlnate  goal  is  to  make  flying  qualities  assessments,  using 
paired  comparisons  between  analog  and  the  two  forms  of  digital  control¬ 
lers*  These  assessments  are  beyond  the  scope  of  the  present  study. 

The  evaluation  of  the  differences  between  baseline  analog  and  digi¬ 
tal  controller  cases  will  constitute  a  three-stage  process*  In  the 
first  stage  the  theoretical  control  law  will  be  tested  against  an  ideal¬ 
ized  analytical  model  of  the  open-loop  YF-16.  In  the  second  stage  the 
control  law  will  be  implemented  in  the  microprocessor-based  system  and 
again  flown  against  an  idealized  model  of  the  open-loop  system.  Here  a 
TR-48  analog  coaqmter  will  serve  as  a  surrogate  YF-16.  In  the  third 
stage  the  hardware-implemented  control  law  will  be  tested  against  a 
simulation  of  the  open- loop  YF-16,  where  now  the  Navion  VRA  will  serve 
as  the  host  aircraft.  This  three-stage  procedure  should  be  very  effec¬ 
tive  in  apportioning  the  errors  introduced  by  the  various  interacting 
elements. 

The  report  is  organized  into  five  sections  and  nine  appendices.  The 
heavy  use  of  appendices  serves  to  increase  the  flew  of  the  development 
by  removing  routine  developments  (e.g.,  first-order  forms  of  the  equa¬ 
tions  of  motion)  and  detailed  theoretical  proofs  from  the  main  body  of 
the  report. 

Section  II  give 8  a  brief  description  of  the  task  and  an  overview  of 
the  ''baseline"  analog  experiment,  the  WLT  direct  side  force  mode.  Sec¬ 
tion  III  treats  the  ZOH  control  law  synthesis  and  validation  in  much  the 
same  manner  as  is  done  for  the  slewer  in  Section  IV.  The  report  con¬ 
cludes  with  Section  V,  conclusions  and  recommendations. 


SECTION  II 


TASK  DESCRIPTION  AND  EQUATIONS  OF  MOTION 


A.  INTRODUCTION 

This  section  focuses  on  a  description  of  the  baseline  analog  wings- 
level  turn  (WLT)  and  those  facets  of  the  testing  approach  that  are  a 
necessary  background  for  understanding  the  digital  implementation. 
There  is  a  brief  description  of  the  baseline  analog  task  and  experiment, 
together  with  a  description  of  the  feedforward  and  feedback  filter  sec¬ 
tions  used  to  cross-couple  the  controller  and  modify  the  open-loop 
dynamics  of  the  simulator  aircraft. 

First-order  forms  for  the  equations  of  motion  are  given  in  the 
s-  and  z-domains  and  for  the  ZOH  and  slewer. 

B.  DESCRIPTION  OF  TASK 

The  purpose  of  the  baseline  analog  flight  test  phase  of  the  Refer¬ 
ence  1  program  was  to  fill  gaps  in  the  data  base  as  required  to  develop 
a  handling  qualities  criterion.  A  tight  tracking  task  was  required  to 
separate  good  and  bad  configurations,  i.e.,  to  force  the  pilot  to  maxi¬ 
mum  effort  and  thereby  expose  deficiencies  which  may  not  otherwise  have 
been  evident. 

The  primary  task  selected  in  Reference  1  was  air-to-air  tracking. 
This  task  was  ideal  because  the  target  motions  could  be  tailored  to 
exercise  a  broad  spectrum  of  frequencies  in  the  tracking  aircraft. 

The  control-configured  vehicle  (CCV)  mode  selected  was  the  wings- 
level  turn  (Ay)  mode.  This  mode  was  selected  because  it  has  consider¬ 
able  potential  for  future  CCV  aircraft  for  air-to-air  and  air-to-ground 
applications.  The  approach  of  Reference  1  was  to  generate  a  series  of 
configurations  with  adverse  and  proverse  roll  and  yaw  coupling  in  the 
wings-level  turn  mode.  These  configurations  were  designed  to  test  the 
bandwidth  hypothesis;  that  is  to  say,  if  the  hypothesis  is  valid  then 
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the  pilot  ratings  should  correlate  with  bandwidth  regardless  of  the  type 
of  coupling*  Based  on  this  line  of  reasoning*  the  following  configura¬ 
tions  were  developed: 

1)  Wings-level  turn  with  adverse  and  proverse  yaw 
coupling  designed  to  vary  heading  bandwidth  from 
zero  to  nearly  7  rad/sec* 

2)  Wings-level  turn  configurations  with  adverse  and 
proverse  roll  coupling*  designed  to  give  the  same 
heading  bandwidth  as  the  configurations  in  Item  1* 

If  the  bandwidth  hypothesis  is  valid*  these  configurations  with  similar 
values  of  heading  bandwidth  should  receive  similar  pilot  ratings  and 
commentary* 

C.  DESCRIPTION  OF  BASELINE  ANALOG  EXPERIMENT 

The  tracking  kinematics  for  the  ideal  case  (no  coupling  or  uncan¬ 
celled  aircraft  modes)  in  the  wings-level  turn  mode  are  summarized  in 
Figure  1*  The  block  diagram  in  Figure  1  indicates  the  interrelationship 
between  the  Idealized  aircraft  dynamics*  the  air-to-air  tracking  kine¬ 
matics*  and  the  target  heading,  The  tracking  kinematics  appear  in 

the  feedback  transfer  functions  of  this  block  diagram  and  result  in  a 
numerator  zero  at  UQ/R  (aircraft  speed/range)*  The  effect  of  this  zero 
on  the  piloted  loop  closure  in  attempting  to  null  the  aim  error,  e,  is 
shown  in  the  root  locus  plot  at  the  bottom  of  Figure  la.  The  closed- 
loop  characteristic  equation  is  seen  to  be  well  damped  when  UQ/R  is 
small,  that  is,  at  large  values  of  range.  At  low  values  of  range*  or 
when  Uq/R  Is  large,  the  closed-loop  characteristic  roots  are  seen  to  be 
lightly  damped.  Physically  this  stems  from  the  fact  that  e  is  primarily 
set  by  heading  when  UQ/R  is  large;  whereas,  when  Uq/R  is  small,  e  is 
strongly  affected  by  lateral  displacement  (which  Involves  an  additional 
integration).  Formation  flying  represents  the  limiting  case  of  this 
with  R  “  0. 

Due  to  structural  limitations  of  the  side  force  generators,  the 
Princeton  University  Variable  Research  Aircraft  (VRA)  has  a  maximum 
maneuvering  speed  of  105  kt  -  well  below  typical  air-to-air  combat 


speeds*  It  was  therefore  necessary  to  adjust  the  range  between  the 
target  and  attacker  in  our  experiment  to  make  the  parmeter  UQ/R  con¬ 
sistent  with  a  typical  air  combat  encounter*  The  effect  of  range  on 
U0/R  for  our  test  conditions  (105  kt  TAS)  is  shown  in  Figure  lb*  Here 
we  can  see  that  typical  combat  parameters  of  M  -  0*86  at  20,000  ft  and  a 
range  of  600-1200  yards  converts  to  100-200  yards  at  the  VRA  testing 
apeed  (105  kt)*  Shorter  ranges  result  in  relatively  large  values  of 
Uq/R,  e*g*,  values  which  approach  the  piloted  crossover  region  in  the 
vicinity  of  1  rad/sec*  The  effect  of  this  on  sight  dynamics  is  shown  in 
Figure  2  for  the  baseline  configuration  (WLT1)  used  in  the  flight  test 
experiment*  The  frequency  response  phase  plot  in  Figure  2  indicates 
that  for  frequencies  below  UQ/R  the  pipper  error  dynamics  are  very 
lightly  damped,  whereas  for  a  large  range  of  frequencies  well  above  u0/R 
the  pipper  error  dynamics  are  equivalent  to  the  heading  dynamics. 
Hence,  for  values  of  range  where  UQ/R  is  well  below  the  piloted  cross¬ 
over  frequency.  It  is  appropriate  to  use  heading  as  the  controlled 
variable  when  applying  the  bandwidth  hypothesis*  Tracking  at  close 
ranges,  where  UQ/R  is  large  enough  to  be  near  the  region  of  piloted 
crossover  (I  rad/sec),  was  found  to  be  impractical  during  initial  flight 
test  evaluation  because  of  the  very  light  damping  of  the  pipper  error 
dynamics*  The  formal  runs  were  conducted  so  that  the  safety  pilot  had 
control  over  range  and  maintained  the  target  aircraft  at  a  nominal  range 
of  150  yards  throughout  the  data  runs*  This  was  accomplished  by  using  a 
series  of  concentric  range  circles  painted  on  the  aircraft  windscreen 
and  sized  so  that  the  target  aircraft's  wingspan  would  be  coincident 
with  the  target  circle  at  a  range  of  150  yards* 

The  primary  disadvantage  of  testing  at  speeds  well  below  M  ■  0.8  is 
that  it  is  not  possible  to  correctly  simulate  the  aircraft  dynamics  and 
the  pilot  acceleration  cues  simultaneously.  This  may  be  seen  from  the 
equation  for  lateral  acceleration: 

®yrt>  *  U0(8  +  r)  -  g* 
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Figure  2.  Comparison  of  Heading  Dynamics  and  Pipper 
Error  Dynamics  for  Wings  Level  Turn 
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Figure  3.  Ideal  Pilot  Input  Signal  to  Roll  Autopilot 

If  the  3  and  r  responses  are  correct,  the  lateral  acceleration  will  be 
scaled  down  by  the  Inertial  speed  UQ.  In  the  present  experiment  we 
elected  to  maintain  the  Integrity  of  the  sideslip  and  yaw  rate  responses 
at  the  expense  of  side  acceleration  cues,  which  were  about  a  factor  of  5 
less  than  those  corresponding  to  M  ■  0*80.  This  was  done  In  accordance 
with  the  notion  that  visual  cues  are  more  dominant  than  acceleration 
cues  In  air-to-air  tracking;  and  with  the  VRA's  maximum  lateral  acceler¬ 
ation  (0*5  g)  capacity*  Lateral  accelerations  as  high  as  0*5  g  were 
utilized  frequently  during  the  experiment*  This  would  translate  to 
about  2.5  g  at  M  -  0. 8. 

The  air-to-air  tracking  scenario  was  developed  to  maximize  the  prob¬ 
ability  of  exposing  deficiencies  in  the  tracking  aircraft.  This  was 
accomplished  by  controlling  the  target  aircraft  heading  0|>2  in  Figure  1) 
in  a  random-appearing  fashion  which  resulted  in  a  power  spectrum  concen¬ 
trated  in,  but  evenly  spaced  over,  the  frequency  range  of  interest* 
I*  A.  M*  Hall  developed  such  a  signal  for  the  purpose  of  identifying  the 
frequency  response  characteristics  of  aircraft  in  flight*  The  signal 
used  is  shown  in  Figure  3*  The  signal  In  Figure  3  was  played  through 
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the  target  aircraft  lateral  autopilot  via  a  switch  controlled  by  the 
target  aircraft  pilot*  This  signal  resulted  in  approximately  three- 
quarters  of  full  aileron  travel  at  the  testing  speed  of  105  kt*  The 
target  aircraft  was  maintained  at  constant  altitude  during  the  run.  It 
was  intended  to  utilize  variations  in  the  input  series,  such  as  playing 
it  backwards  or  from  the  middle  to  the  ends,  etc.;  however,  the  evalua¬ 
tion  pilots  felt  that  the  task  remained  unlearned  (random)  and  therefore 
such  variations  were  not  utilized. 

D.  BASELINE  ANALOG  CONFIGURATION 

The  Princeton  University  Variable  Stability  Research  Aircraft  (VRA) 
is  a  fly-by-wire  response  feedback  simulation  utilizing  hydraulically 
actuated  controls.  These  controls  include  flaps  which  move  up  as  well 
as  down  and  side  force  generators  as  sketched  in  Figure  4.  A  block  dia¬ 
gram  of  the  VRA  as  mechanized  for  the  in-flight  simulation  in  this 
program  is  shown  in  Figure  5.  The  Cfi  and  CF  matrices  in  Figure  5  were 
calculated  to  allow  the  VRA  to  respond  like  the  YF-16  at  a  flight  condi¬ 
tion  of  M  *  0.8  at  20,000  ft.  The  generic  variation  of  roll  and  yaw 

coupling  in  the  flight  test  experiment  was  achieved  via  the  aileron  and 
rudder  crossfeed  boxes  in  Figure  5. 

E.  IN-FLIGHT  VEHICLE  IDENTIFICATION 

A  primary  problem  with  much  of  the  DFC  data  generated  to  date  is 
that  the  actual  controlled  element  tested  was  not  quantitatively 
defined.  To  avoid  any  uncertainties  in  defining  the  controlled  element 
for  each  configuration  tested  in  this  experiment  it  was  decided  to  run  a 
frequency  sweep  between  the  DFC  input  and  the  aircraft  heading  which 
could  then  be  Fast  Fourier  Transformed  to  obtain  the  frequency  response 
directly.  This  technique  had  the  secondary  advantage  of  determining 
whether  or  not  it  is  practical  to  write  flying  qualities  criteria  in 
terms  of  frequency  response  characteristics.  The  method  for  generating 
the  frequency  sweep  was  extremely  simple  in  that  the  pilot  simply  exer¬ 
cised  the  DFC  control  (rudder  pedals)  at  ever-increasing  frequency  dur¬ 
ing  a  single  run.  Rudder  pedal  input  and  output  yaw  rate  were  recorded 
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and  Fast  Fourier  Transformed  with  excellent  results,  i.e.,  very  little 
data  scatter  in  the  frequency  range  of  interest.  An  example  of  a  Bode 
plot  generated  in  this  manner  is  shown  in  Figure  6. 

F.  VBA  SETUP  FOB  ANALOG  DIRECT  FORCE 

TESTS  OF  REFERENCE  1 

The  VRA  setup  for  the  analog  simulation  of  the  YF-16  wings-level 
turn  is  shown  in  Figure  7.  In  Figure  7,  S  refers  to  the  simulator 
(Navion)  while  M  refers  to  the  model  (YF-16). 

The  Cp  and  Cg  matrices  are  a  hardwired  part  of  the  analog  control 
system;  there  is  a  pot  to  set  for  each  element.  However,  the  matrices 
do  contain  zero  entries  (e.g.,  $  +  Sr,  p  ♦  Sgp).  Some  of  these  were  not 
critical  in  the  analog  flight  tests  and  were  left  at  zero. 

The  "mixing  box"  and  Cp  matrix  in  Figure  7  need  some  explanation, 
which  requires,  as  a  prerequisite,  a  look  at  the  lateral  VRA  system 
prior  to  the  Reference  1  study  (refer  to  Figure  8) .  In  Figure  8  the 


w (rod /sec) 


Figure  6.  Fourier  Transformed  Heading  Responses, 
Configuration  WLT1,  Minimal  Coupling 
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6S  -  («r  (rudder),  6SF  (side  force),  <SA  (aileron) )Navlon 

6SF* 

■  Pilot  input  (rudder  pedal  deflection  in  Navion) 

Cg  “  Response  feedback  (gain)  matrix 

Cp  -  Feedforward  gain  matrix  (control  crossfeed) 

Figure  7.  Analog  Test  Setup 


This  matrix  has 
too  aany  zeros 


Figure  8.  Lateral  VRA  System,  Prior  to  Ref*  1  Study 
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limited,  existing  control  coupling  gain  matrix  was  not  utilized.  The 
feedforward  gain  matrix  had  too  many  zero  entries  to  be  useful;  to 
remedy  this  the  VRA  personnel  Installed  a  "mixing  box"  to  furnish  the 
necessary  feedforward  gains.  This  box  furnished  a  matrix  of  gains  (a 
1  x  1  matrix)  which  interfaced  with  a  1  *  3  matrix  ,tbox"  furnished  by 
STI.  This  box,  in  turn,  coupled  with  another  3*3  matrix  of  gains,  Cp, 
furnished  by  Princeton.  The  scalar  was  utilized  as  a  pure  gain  dur¬ 
ing  the  course  of  the  study;  although  provision  was  made  for  switching 
in  an  integrator  (Km  *  K  /s). 

The  purpose  of  the  STI-fumished  box  (Cp)  was  to  provide  control 
crossfeeds  that  were  functions  of  frequency;  to  be  specific,  a  1  *  3 
matrix  of  first-order  filter  sections. 

With  this  configuration,  the  necessary  control  crossfeeds  for  the 
various  YF-16  CCV  diodes  could  be  implemented  as  filter  sections  (rather 
than  Just  pure  gains)  and  thus  satisfy  the  frequency-dependent  coupling 
numerator  requirements  of  Reference  1. 

For  the  present  study,  two  possibilities  exist  as  far  as  the  feed¬ 
forward  matrix  is  concerned: 

1)  Princeton  can  expand  the  "old"  CF  matrix  (with 
the  many  zero  entries)  to  a  more  usable  form. 

2)  The  mixing  box  can  be  reinstalled  in  the  VRA  to 
give  a  complete  gain  cross-coupling  (feedforward) 
capability. 

Approach  2  was  the  one  ultimately  implemented. 

The  filter  break  frequencies  are  very  close  for  the  wings-level  turn 
and  hence  the  crossfeeds  can  be  approximated  with  pure  gains. 

Another  important  point  concerning  the  Reference  1  test  is  that  the 
analog  gains  (furnished  by  Cg  and  the  "mixing  box"/STI  Cp  filter  matrix) 
were  utilized  to  force  the  bare  airframe  Navion  to  give  an  equivalent 
lower-order  model  of  a  compensated  YF-16  configuration  (that  is,  a  bare 
airframe  YF-16  with  all  the  relevant  feedforward/feedback  loops  closed). 
The  term  "lower-order  model"  is  used  because  the  YF-16,  with  compensa¬ 
tion  networks,  is  a  relatively  high-order  system.  The  lower-order 


equivalent  model  was  fourth  order,  which  leads  to  a  considerable  simpli¬ 
fication  in  the  design  of  the  wlngs-level  task  (see  Section  H). 

6.  TP-16  BAKE  AHFBAMB  CONSIDERATIONS 

Using  the  VRA's  analog  fly-by-wire  capability  to  go  directly  from 
the  bare  airframe  Navion  to  a  compensated  YF-16  is  a  reasonable  approach 
for  an  analog  fly-by-wire  mechanization.  For  the  digital  tests  this 
process  would  obscure  the  artifacts  Introduced  by  the  zero-order  hold, 
variation  in  the  data  rate,  etc.  For  the  digital  tests  it  is  desirable 
to  utilize  the  analog  fly-by-wire  system  to  force  the  Navion  to  appear 
as  a  bare  airframe  YF-16.  Then  the  Micro-DFCS  can  be  used  to  implement 
the  digital  controllers.  Thus  the  implementation  of  Figure  9,  where  Cp 
and  Cg  are  "analog"  gain  matrices  (furnished  by  some  combination  of 
either  existing  VRA  capabilities  or  in  combination  with  the  mixing  box) 
will  allow  an  exact  replication  of  the  bare  airframe  YF-16  by  the 
Navion.  The  context  of  "exact"  must  wait  upon  the  expository  material 
of  a  later  section. 


Bare  Airframe  YF-16 


Figure  9.  "Analog"  Test  Configuration 
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H.  FIRST-ORDER  FORM  OF  EQUATIONS  OF  MDTIOR 


Let  the  time  invariant  model  of  the  aircraft  to  be  simulated  (YF-16) 
be 

XM  “  fmxm  +  CD 

while  the  simulation  aircraft  dynamics  are  modeled  as 

*S  ”  FSXS  +  GS6S  (2) 

•  e 

The  components  of  x  are  yaw  rate  sideslip  (3),  roll  rate  (<{>),  and 

roll  angle*  The  control  components  are  rudder  (6  ),  side  force  (6  ) 

and  aileron  (6^) • 

Given  the  stability  axis  state  vector,  the  F  matrix  has  the  generic 


NB  Np  0 


<Yr/V)-l  Yb/V  Yp/V  g/V 


Lb  Lp  0 


For  the  three  lateral  controllers,  the  control  effectiveness  matrix 
has  the  generic  form 


8«r  n«sf  n«a 


y5r/v  y6sf/v  ySa/v 


L$r  L6sf  Lj a 


I 

! 


P  la  a  4  x  4,  but  the  control  effectiveness  matrix  Is  a  4  x  3  (not 
square).  The  entries  of  Equations  3  and  4  are  tabulated  in  Appendix  A. 

I.  DISCRETE  FIRST-ORDER  FORM  (ZOH) 

The  discrete  version  of  Equation  1  (or  Equation  2)  has  the  format 

xk+l  "  *xk  +  r6k  (5) 

where 

*  -  £-1[ is  -  p]“  (6) 

t-T 

r  ■  l  -1  f  *>'L)c 

♦(t)  and  T (T)  are  tabulated  In  Appendix  A. 

J.  DISCRETE  FIRST-ORDER  FORM  (SLEMER) 

The  discrete  version  of  Equation  1  Is 

xk+l  *  ♦<T>xk  +  I* 25k-l  (8) 


where 


♦  *  jCl  “  F]”1 


(9) 


(10) 


r2  "  rZ0H  -  rl  (See  Eq.  7  for  r Z0H ) 


(ID 


SECTION  III 


ZOH  DIGITAL  CONTROLLER 


A.  INTRODUCTION 

A  digital  control  law,  using  a  zero-order  hold  as  the  coupler 
between  the  digital  computer  and  the  control  actuators  is  synthesized 
directly  in  the  z-domaln.  That  is,  the  "emulation"  approach  of  digitiz¬ 
ing  an  existing  analog  control  law  is  avoided.  The  approach  is  to 

develop  discrete  models  of  the  open-loop  YF-16  and  then  apply  an  exten¬ 
sion  of  the  equivalent  stability  derivative  (ESD)  approach  to  synthesize 
the  control  law.  An  overview  of  the  ESD  method,  and  the  extension 
(which  uses  a  pseudo  inverse),  is  described  first.  After  this,  the 
computational  details  of  implementing  the  approach  and  validating  the 
control  laws  are  taken  up. 

An  important  point  concerning  the  analog  baseline  experiment  (Refer¬ 
ence  1)  is  that  analog  gains  were  utilized  to  force  the  bare  airframe 
Navlon  to  give  an  equivalent  lower-order  model  of  a  compensated  YF-16 
configuration  (that  is,  a  bare  airframe  YF-16  with  all  the  relevant 
feedforward/feedback  loops  closed).  The  term  "lower-order  model"  is 
used  because  the  YF-16,  with  compensation  networks,  is  a  relatively 
high-order  system. 

The  VRA's  analog  fly-by-wire  system  was  used  to  go  directly  from  the 
bare  airframe  Navlon  to  a  compensated  YF-16  in  the  analog  tests.  For 
the  digital  tests  it  is  desirable  to  utilize  the  analog  fly-by-wire 
system  to  force  the  Navlon  to  appear  as  a  bare  airframe  YF-16.  Then  the 
Mlcro-DFCS  can  be  used  to  Implement  the  digital  controllers. 

The  point  is  a  crucial  one  and  deserves  elaboration.  If  the  digital 
control  law  is  designed  to  force  the  open-loop  Navion  to  look  like  a 
closed-loop  YF-16  (at  the  sampling  instants),  then  (in  the  inter-sample 
period)  the  Navion  will  respond  like  a  Navion  —  not  the  open-loop 
YF-16.  Consequently,  one  must  first  "wrap"  an  analog  fly-by-wire  system 


around  the  Navion  to  insure  that  the  inter-sample  aircraft  performance 
resembles  that  of  an  open-loop  YF-16.  This  is  especially  critical  if 
one  intends  to  evaluate  the  effect  of  a  relatively  low  data-rate  digital 
control  working  against  a  relatively  wideband  open-loop  system. 

In  view  of  this,  there  is  an  additional  requirement  to  exercise  the 
ESD  method  and  define  the  (continuous)  feedback  and  feedforward  gains 
that  force  the  Navion  to  emulate  the  open-loop  YF-16  model. 

The  validation  procedure  is  a  three-step  process.  First,  the  theo¬ 
retically  derived  control  law  is  tested  against  a  theoretical  continuous 
model  of  the  open-loop  YF-16.  This  was  done  using  the  TOTAL  program 
(Reference  5).  Various  metrics,  such  as  w'-  and  s-plane  comparisons, 
and  frequency  response  and  step  response  comparisons  were  used  to  assess 
the  verisimilitude  of  the  results. 

Second,  the  control  law  was  implemented  into  microprocessor  hardware 
and  validated  against  a  theoretical  model  of  the  YF-16,  that  is,  flown 
against  an  analog  computer  model  of  the  YF-16. 

The  last  step  is  the  actual  flight  test,  with  the  Navion  VRA  simu¬ 
lating  the  open-loop  YF-16.  As  will  be  seen,  the  pseudo  inverse  gener¬ 
alization  of  the  ESD  approach  does  an  excellent  job  of  matching  the 
digitally  controlled  open-loop  YF-16  to  the  low-order  closed-loop  model. 

B.  REVIEW  OF  THE  ESO  METHOD 

An  important  design  approach  to  in-flight  simulation  is  the  equiva¬ 
lent  stability  derivative  (ESD)  technique.  This  approach  can  force  the 
responses  of  a  simulator  aircraft  to  "perfectly"  match  those  of  a  model 
aircraft  if  certain  constraints  pertaining  to  the  number  of  controllers 
and  states  are  satisfied. 

Let  the  time  invariant  (lateral-directional)  model  of  the  aircraft 
to  be  simulated  be 

*m  “  fM%  +  %6M  (I4) 
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while  the  simulation  aircraft  dynamics  are  modeled  as 


xs  -  Fsxs  +  Gs$s  <15) 

The  components  of  x  are  yaw  rate  ($),  sideslip  (8),  roll  rate  ($),  and 
roll  angle  (<j>).  The  control  components  are  rudder  (<5R),  side  force 
(6Sp),  and  aileron  (6^) . 

Given  the  stability  axis  state  vector,  with  F  and  G  matrices  having 
the  literal  form 


Nr 

Nb 

Np 

0 

(Yr/V-1) 

YB/V 

Yp/V 

g/v 

Lr 

lB 

Lp 

0 

0 

1 

0 

n«sf 

N6A“ 

y«r/b 

y«sf/v 

y«a/v 

L6r 

l«sr 

l«a 

0 

0 

0 

(16a) 


(16b) 


it  is  apparent  that  three  controllers  can  provide  "perfect"  model  fol¬ 
lowing  (subject  to  control  actuator  limitations)  since  the  fourth  row  of 
Equation  16  is  an  identity  (it  requires  no  "matching")  for  both  model 
and  simulator* 


Equation  26  can  no  longer  be  forced  to  look  exactly  like  Equation  23 
because  does  not  exist  (it  is  a  4  x  3).  Furthermore,  the  fourth  row 
is  non-zero,  so  that  working  with  a  row-reduced  version  of  T  is  not 
feasible.  However,  the  pseudo  inverse  continues  to  work  and  one  may 
define 


CF  „  (IW1  T'rm  ,  CB  -  (rsrs)-1  rs^m  -  4>s)  (27) 


as  the  best  mean  square  fit. 

To  demonstrate  these  points,  consider  a  two-state,  single-controller 
example: 
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A  straightforward  computation  yields 


Clearly, 


CF 


GsGm  “  “1 

(G'Gg)"1  Gg ($m  -  *8)  -  [-10  -4] 
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and  the  simulator  exactly  matches  the  model 


However,  the  discrete  models  no  longer  exhibit  the  8-plane  identity 
in  the  first  component: 


*«k+l  " 


2e-T  -  e“2T  ]  e-T  _  e-2T 


•2e“^  +  2e”2T  t  -e-"^  +  2e“2T 


1/2  -  e_T  +  (l/2)e“2T 


e-T  _  e-2T 


4e~3T  _  3e“^T 


»-3T  _  „-4T 


— 12e“3T  +  12e“4^  |  — 3e“^T  +  4e“^ 


-1/12  +  (l/3)e“3T  -  (l/4)e“4T 


_e-3T  +  e-4T 


The  pseudo  Inverse  still  yields  an  answer,  but  it  lacks  the  power  o£ 
the  s-plane  result,  since  the  fit  is  no  longer  exact  (it  is  a  mean- 
square  fit).  Thus, 

^au  "  ^F^mu  +  ^Bxsv 


«W‘  rIrB 


<W‘  r;<*„ 


will  not  give  Equation  32  "exactly." 
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D.  A  CHECK  COMPUTATION:  NATION 


CLOSED -LOOP  YF-16 


The  "model"  for  the  digital  in-flight  simulation  is  the  YF-16. 
A  first  step  is  to  wrap  analog  loops  around  the  simulator  aircraft 
(Navion)  so  that  the  digital  control  loops  presume  an  open-loop  YF-16. 
To  demonstrate  the  mathematical  process  we  will  first  present  a  check 
computation  to  verify  the  Reference  1  gains  used  to  force  the  Navion  to 
look  like  a  closed-loop  YF-16.  The  matrices  needed  (see  Appendix  A)  for 
both  the  Navion  and  YF-16  are  tabulated  below.  The  YF-16  flight  condi¬ 
tion  is  for  20,000  ft  and  Mach  0.8. 


Navion 


0.777 

4.68 

-0.432 

0 

•1.0 

-0.3556 

0 

0.172 

1.27 

-12.8 

1 

0 

0 

0 

1 

0 

(34a) 


-6.1  2.41  0.314 

0.0725  0.237  0 

0.77  0  21 

0  0  0 


(34b) 


The  next  set  gives  the  F  and  G  matrices  for  the  YF-16,  with  a  verti¬ 
cal  canard  (CCV  vehicle)  as  per  the  flight  tests  of  Reference  1  (i.e., 
the  analog  flight  tests).  The  3,4  entry  of  the  F  matrix  changes  drama¬ 
tically  from  a  value  of  -0.2077  with  no  canard  to  a  value  of  -25.0  with 
the  canard  (this  is  the  fourth-order  equivalent  model  used  in  the  Refer¬ 
ence  1  tests). 


For  all  cases,  the  state  vector  la 


AX  -  [Ar  AB  AP  A*]' 


(39) 


and  the  control  vector  la 


AS  ■  [ ASj^  ASgp  ASa] 


(40) 


First,  compute  the  gains  (which  determine  the  pot  settings  for  the 
analog  fly-by-wire  computer)  which  force  the  Navion  to  have  the  ESD 
(equivalent  stability  derivative)  model  of  the  YF-16  CCV,  canard,  loops 
closed*  Since  this  is  an  s-plane  computation  we  elect  to  use  Equation 
20*  The  first  part  of  Equation  20,  which  has  the  form 


CB  "  G-l(Fm  -  Fs)  (41) 

is  rewritten  in  a  form  suitable  for  solution  using  either  Cramer's  rule 
or  Gaussian  elimination.  That  is, 

gbcB  “  (Fm  “  Fs) 

This  gives 


-6*  1 

2.41 

1 

CO 

• 

o 

1 

-2.328 

4.237 

0.16 

0 

0.0725 

0.237 

0 

cB  - 

0.0204 

0.0591 

0 

-0.13405 

Gaussian  elimination  then  gives  (nine-place  accuracy  is  used  only  to 
facilitate  checking  the  answer): 


”  0.359857139 

-0.456220933 

-0.019653205 

-0.144939218" 

-0.024006931 

0.388928345 

0.006012056 

-0.521273868 

(43) 

_  0.238900576 

-1.645652851 

-0.081612716 

-1.  185161753  _ 

confute  Cp,  the 

second  part  of 

Equation  20 

is  also  written 

In  a 

form  suitable  for  Gaussian  elimination. 


G8CF 


G 


m 


"-6.1 

2.41 

-0.314" 

’-3.925 

4.766 

1.897  ' 

0.0725 

0.237 

0 

Cp  “ 

0.02988 

0.01802 

-0.0336 

0.77 

0 

21  . 

7.178 

5.655 

49.6 

The  solution  is 


CF 


0.603820992 

-0.058637225 


-0.683783153 

0.285207927 


-0.436629058 

-0.008204191 


(45) 


0.319669421 


0.294357763 


2.37791449 


These  gains  agree  with  the  Reference  1  results.  Next,  compute  the 
gains  for  Navion  to  bare  airframe  YF-16.  Since  the  previous  sample  cal¬ 
culations  document  the  solution  procedure  for  solving  Cg  and  C^,  we  will 
dispense  with  the  details  and  present  the  results. 


”-0.051566838 

-0.653322784 

-0.066549239 

-0.198437737” 

0.01577466 

0.409138817 

0.020357890 

-0.501321790 

(46) 

_  0.020938403 

-1.747473404 

0.205773472 

0.  007276050_ 

0.641584848 

-0.661386936 

-0.478658047 

-0.062088192 

0.282913725 

-0.004207981 

(47) 

_  0.339808556 

0.306631807 

2.355646033 

Next,  compute  the  digital  gains  which  force  the  bare  airframe  to 
match  the  closed-loop  YF-16  (at  the  sampling  instants).  First,  the 
plant  equations  must  be  properly  discretized,  taking  into  account  both 
the  sample  rate  and  the  data  hold.  This  is  done  in  the  next  section, 
where  the  focus  is  on  the  ZOH. 

E.  ZOH  ESD  CONTROL  LANS 

To  review,  the  Navion' s  analog  system  will  be  utilized  to  force  the 
Navion  to  appear  as  the  open-loop  (bare  airframe)  YF-16  with  canard. 
(In  Figure  11,  the  subscript  C  refers  to  the  continuous  gain.)  In  turn, 
the  Navion's  Micro-DFCS  is  then  used  to  implement  feedforward  and  feed¬ 
back  gains  to  make  the  open-loop  YF-16  appear  as  a  closed-loop  YF-16 
(Figure  12) • 

The  task  now  is  to  find  the  ESD  gains  given  finite  T.  Recall  that 
the  continuous  gains  were  found  very  easily;  given  three  controllers  and 
four  states  (one  being  an  identity),  exact  matching  was  possible.  For 
finite  T,  this  fortunate  situation  no  longer  pertains;  the  fourth  equa¬ 
tion  is  no  longer  an  identity.  This  is  best  illustrated  with  the  aid  of 
Table  1,  which  describes  the  ZOH  case  (T  -  0.1,  YF-16  CCV  canard,  bare 
airframe) : 
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'  \ 
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1 

TABLE  1.  s-  AND 

z-PLANE  OPEN  LOOP 

PLANT  REPRESENTATIONS 

•1 

i 

x  “  Ax  +  B6 

2 

|  AHAT(ROU 

1  >• 

-.4318900000 

10.20000000 

-.4i6eee00eeE-0i 

0. 

1  AflAT  ( ROU 

2  >• 

-1.000000000 

-.3060000000 

e. 

.38800000006-01 

N  AflAT  (ROW 

3)- 

1.670000000 

-50.00000000 

-2.330000000 

N  APIAT(ROU 

4>- 

0. 

0. 

l  .000000000 

0. 

ft  Tiwmrair 

TT- 

-4. i7oeeo000 

4.62000000® 

"2'.  170'000W0~ 

Ft  BHAT  ( ROU 

2)- 

.31 8O000O00E-01 

.l9ioeeee00E-0t 

-.3570ee0OO0E-01 

ft  BP1ATCR0U 

3)- 

7.630000000 

5.930000000 

49.10000000 

P  BflAT  ( ROU 

4»- 

0. 

0. 

0. 

1 

xk+l  -  Fxk  +  G6k 

S  mart  rou 

1  )- 

. 90856 IT123 

•975S276552 

- . 3S02672856E-02 

. 19267557606-02 

S'  rnATtRou 

2  >• 

-.94716629936-01 

.9204149530 

.36268690486-03 

.37566966706-02 

S  FUAT (ROU 

3)- 

.3670070570 

-4.234745790 

.7912792594 

-.87163153416-02 

S  FT1AT  ( ROU 

8  CHAT (ROU 

4  )* 

. 1524291762E-01 

-.2246473026 

.8917931248E-01 

.9997011628 

-74010594654' 

.'4443513625 

71978434172 

1  CHAT (ROU 

2)- 

.2334360234E-01 

-. 20423525 3 IE-01 

-. 133427171  IE-01 

l  CHAT (ROU 

3  )• 

.6097314035 

.5949648390 

4.419801283 

»  CHAT (ROU 

h 

C 

4>- 

<  3320928555E-0 1 

• 29422492SSE-01 

.2287260839 

- 

the  identities  in  the  fourth  row  of  the  A  and  B  matrices  (0,  0, 

l.  o 

and 

0,  0,  0).  The 

fourth  row  of  F  and 

G  must  take  on 

the  numerical 

r 

f  values  dictated  by 

1 

♦  (T)  -  *C_1  [is  - 

-1 

A] 

(48) 

i 

1 

t-T 

n  and 

T 

r<T)  -  f  <Kt  - 
•/o 

S)B  d5 

(49) 

j  and 

hence  the  identity 

is  lost*  Thus  the  time  domain 

procedure  of 

1  throwing 

away  the  fourth  equation,  to  obtain  a  reduced  control  effec- 

jj  tiveness 

matrix  that  is 

invertible,  is  no  ! 

Longer  possible. 

To  preserve 

|  the 

ESD 

concept  we  may  appeal  to  the  pseudo 

-inverse,  and  then  check,  via 

I 

to  see  how  good  the  result  is*  The 

pseudo-inverse 

equations  are 

1  described  next* 
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^  V  *  ,  -  \  - 

Let  the  time  domain  equations  for  the  open-loop  YF-16  (the  "simula¬ 
tor/1  subscript  8)  be 

xa  -  (IMAT)sxs  +  (JMAT)a6s  OL  YF-16  (50) 

Let  the  time  domain  equations  for  the  closed-loop  YF-16  (the  "model," 
subscript  m)  be 

^  “  (XMADaXn,  +  (YMAT)n6m  CL  YF-16  (51) 

Use  Option  87  to  get  the  discrete  (ZOH)  form: 

Xfc+i  -  (LMAT)xk  +  (MMAT)6k  CL  YF-16  (Model)  (52) 

x^i  -  (NMAT)xk  +  (OMAT)6k  OL  YF-16  (Simulator)  (53) 

Since  OMAT  is  the  control  effectiveness  matrix  of  the  "simulator"  (OL 
YF-16),  we  need  its  pseudo-inverse.  Let 

PHAT  -  (OMAT)t  -  [O'  0]_1  0'  (54) 

where  t  indicates  the  pseudo-inverse. 

Clearly,  the  feedforward  gain  matrix,  Cp,  is  given  at 

QMAT  -  CF  -  (PMAT)*(MMAT)  (55) 

In  a  like  manner,  the  feedback  gain  matrix,  Cfi,  is  given  by 

UMAT  -  CB  -  (PMAT)ILMAT  -  NMAT]  (56) 


*For  a  description  of  TOTAL  matrix  notation  and  analysis  options, 

see  Reference  5. 
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Equations  55  and  56  define  the  gain  matrices.  We  oust  now  check  by  sub¬ 
stituting  Equations  55  and  56  back  into  the  equations  for  the  open-loop 
system  (Equation  53)  and  see  how  well  the  result  compares  with  the  model 
(Equation  52). 

G.  COMPARISON  EQUATIONS  —  STATE  TRANSITION 

With  the  control  law 

«k8  “  CBxkg  +  =  (UMAT)xkg  +  (QMAT)6kin  (57) 

Equation  53  becomes 

xk+1  -  [NMAT  +  OMAT*UMAT]xs  +  (OMAT)(QMAT) 6m 

-  RMATxs  +  SMAT6m  (58) 

Since  the  model  is 

Xkm  -  (LMAT)Xkm  +  (MMAT)6km  (59) 

compare  RMAT  with  LMAT  and  SMAT  with  MMAT  (this  is  done  in  Appendix  C). 

H.  w'  —COMPARISON 

Given 

xk+l  M  *f>xk  +  r6k 
take  the  z-transform  to  form 

[Iz  -  $]x(z)  =*  T6k  +  zx(0) 


(60) 


(61) 
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Substituting  and  clearing  gives  [for  convenience,  drop  the  prime  on  w' 
and  set  x(0)  -  0] , 

Iv  +  ■£  (I  +  +>“1  (1  -  *)]x(v)  -  *  [I  +  2^  +  l)r6(v)  (62) 

First,  generate  the  v'  format  for 

Xfc+i  -  RMATxk  +  SMAT6m 

We  can  define  this  to  be 

|lw  +  y  (I  +  RMAT)-1  (I  -  RMAT )  X  -  |  t1  “  RMAT]“1SMAT  +  1^  6  (w) 

(63) 

[Iw'  +  WMATJX  -  EMAT6m  j/T  +  l)  (64) 

This  can  be  compared  with  the  continuous  model.  That  is,  compare 

(Is  -  XMAT)X(s)  -  [YMAT] 6S  (65) 

against 

(Iw  +  WMAT)X(w)  -  [EMAT] 6a  (66) 

See  Appendix  C  for  a  specific  comparison. 

I.  MACROS 

Three  macros  were  used  to  implement  the  previously  described  compu¬ 
tational  algorithms.  The  first,  AKEY,  sets  up  the  state  transition 
model,  from  the  time  domain  model,  for  both  the  simulator  (Equation  50) 
and  the  model  (Equation  51). 


OPTION  >CREATE y rtKEY y COPY  y XMAT  y AMAT y COPY  y  YMAT y BMAT  y 87  y COPY  y PH A f y LMAT t 


>COPY * GMAT y MMAT » COPY y IMAT y AMAT » COPY y  JMAT y BMAT y 87 y 
>COPY y PMAT y  NMAT y COPY  y GMAT y OMAT 

Next,  BKEY  computes  the  pseudo-inverse  and  the  feedback/feedforward 
gains.  It  lists  LMAT,  MMAT,  NMAT,  OMAT,  PMAT,  QMAT,  UMAT. 

OPTION  > CREATE y  BKEY r COPY  > OMAT y AMAT  r 76 y COPY y CM AT r PMAT y  COPY » CM AT y  AMAT» 

>COPY f  OMAT y  BMAT  y  7 A  *  COPY y  CMAT * AMAT  t 75  y  COPY » CMAT y  AMAT  t COPY  t PMAT  y  BMAT  y 

>74 y  COPY  y  CMAT y  PMAT  y  COPY  y  PMAT  y  AMAT  y  COPY y  MMAT y  BMAT  y  74  y  COPY  y  CMAT  y QMAT  y 

>COPY y  LMAT y  AMAT  y  COPY  y  NMAT  y  BMAT  y  73  y COPY  y  CMAT y BMAT  y 

>COPY y  PMAT  y  AMAT  y 74  y  COPY  y  CMAT  y  UMAT  y  LMAT  y  MMAT  y  NMAT  y OMAT y  PMAT  y 

>QMAT y  UMAT 

Next  use  CKEY  to  check  the  results.  Cg  and  Cy  are  substituted  into 
the  "simulator"  equations,  and  this  result  is  then  compared  against  the 
"model"  equations. 

OPTION  >CREATE t CKEY  y  COPY  y  OMAT  y AMAT  y  COPY y  UMAT  y  BMAT  y  74  y  COPY y CMAT  y  BMAT  y 

>COPY y  NMAT  y AMAT  y  72  y  COPY  y  CMAT  y  RMAT y  COPY  y  QMAT  y  BMAT  y  COPY y  OMAT  y  AMAT  y 

>74  y  COPY  y  CMAT y  SMAT  y  COPY  y  RMAT  y  BMAT y  COPY y  UMAT y AMAT  1 72 r COPY y CMAT  y  AMAT  y 

>75  y  COPY  y  CMAT  y  BMAT  y  COPY  y ZMAT y  AMAT  y  74  y COPY  y EMAT  y  TMAT  y COPY y  UMAT  y  AMAT  y 

>COPY y  RMAT  y  BMAT  y  73  y  COPM - CMAT  y  BMAT  y  COPY y TMAT  y AMAT y  74  y  COPY  y  CMAT r  UMAT y 

>COPY y  TMAT  y  AMAT  y  COPY  y  SMAT  y  BMAT  y  74  y  COPY  y  CMAT  y  EMAT  y  COPY  y  WMAT y  BMAT  y 

>COPY y  HMAT y  AMAT  y  74  y  COPY y  CMAT  y  AMAT  y  7 1 

An  annotated  run  for  the  10  Hz  case  is  given  in  Appendix  C. 
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The  Cjj  and  C?  gains  are  tabulated  in  Table  2.  The  printout  rounds 
the  gains  to  four  figures,  although  the  Internal  precision  is  much 
higher*  To  check  the  sensitivity  to  roundoff,  we  reran  CKEY  for  the 
5  Hs;  and  10  Hz  cases  using  the  four-place  numbers  printed  in  Table  2* 
There  were  no  significant  changes  in  the  closed-loop  plant  description, 
indicating  four-place  accuracy  is  sufficient* 

K.  VALIDATION:  THEORETICAL  CONTROLLER, 

TOTAL  MODEL  OF  YF-16 

The  first  method  for  validating  the  theoretical  controller  against 
the  theoretical  YF-16  model  was  outlined  in  Section  H.  For  the  10  Hz 
case,  CB  and  Cp  are  substituted  into  the  "simulator"  equations,  result¬ 
ing  in  the  closed-loop  format 

(Is  -  XMAT)X(s)  -  [YMAT]5S  (67) 

This  can  be  compared  against  the  YF-16  s-plane  reference  model,  given  by 

(Is  +  WMAT)(X(w)  -  [EMAT]6S  (68) 

This  is  carried  out  in  Table  3*  As  can  be  seen,  even  at  10  Hz  the 
pseudo  Inverse  has  done  an  excellent  job  of  matching  the  digitally  con¬ 
trolled  open-loop  YF-16  to  the  desired  closed-loop  YF-16  configuration. 

The  second  method  of  validation  is  to  compare  the  frequency  response 
of  the  continuous  variables  of  the  discretely  controlled  continuous 
system  against  the  frequency  response  of  the  baseline  closed-loop  YF-16. 
A  typical  (10  Hz)  set  is  given  in  Figures  13  and  14.  The  matches  are 
extremely  good  at  low  frequencies  and  exhibit  the  properties  one  expects 
of  digitally  controlled  systems  at  higher  frequencies.  A  more  compre¬ 
hensive  set  (10  Hz  and  50  Hz)  is  listed  in  Appendix  D.  The  third  method 
for  comparison  is  to  match  step  responses  (a  step  of  pedal  input) . 
These  are  extremely  close  and  a  comparison  of  printed  data  is  required 
in  order  to  see  the  differences  (see  Table  4). 


TABLE 


w 

CO  cm 

0>  co 

8  CM  O 
Os 


3780E-02  -.2687E-02  i  .9752  -.3721E-02  .2499E-01  !  — 1009  -.4136 

- - - j - 1 - 

.9388  -. 3329E-01  .6451E-01  .6400  .3067  -.7433E-01  -.4046 
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'59301931  0.000040843  0.000011374  0.000006663  0.000001716-0.000500803-0.000431708 
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L.  VALIDATION:  MICROPROCESSOR-BASED  CONTROL  LAN, 
ANALOG  COMPUTER  MODEL  OF  YF-16 


The  ZOH  control  law  was  Implemented  in  the  microprocessor  and  val¬ 
idated  on  the  analog  computer  using  an  "ideal"  aircraft  model.  The 
wings-level  turn  was  successfully  flown.  A  comparison  between  the  theo¬ 
retical  response  and  digital  implementation  in  the  hybrid  simulation  is 
given  in  Figure  15. 

As  can  be  seen  by  inspection,  the  response  characteristics  of  the 
micro-processor  coupled  controller  closely  match  those  of  the  continuous 
model.  The  magnitude  differences  between  the  two  sets  of  curves  were 
due  to  scaling  limitations  imposed  by  the  coupling  betwen  the  micro¬ 
processor  and  the  TR-48  bare  airframe  simulation.  The  more  favorable 
scaling  dynamics  of  the  coupling  between  the  microprocessor  and  the 
Navion's  analog  system  assured  us  that  the  ZOH  would  function  appropri¬ 
ately  under  actual  flight  test. 


M.  VALIDATION:  MICROPROCESSOR-BASED  CONTROL  LAN, 
NAVION  AS  YF-16  SIMULATOR 


The  microprocessor-based  control  law  was  Installed  in  the  Navion 
Variable  Stability  Aircraft  and  a  wings-level  turn  successfully  flown 
(see  Figure  16).  The  pedal  was  ramped  in  and  held  at  about  70  percent 
full  deflection.  A  wings  level,  essentially  zero  sideslip  turn  near  the 
Navion  limit  of  ay  i  1/2  g  was  observed  in  flight.  That  the  turn  is 
wings  level  is  evidenced  by  the  roll  rate  trace  in  Figure  16  which 
remains  near  zero  throughout  the  maneuver.  For  a  steady  1/2  g  turn 


.  0.5g  0.5*32.2  ft/sec^  *  c,  _  ,  ,  , 

r  ■  "  ■  —177  ft/sec  *  57 -3  de*/rad 


*>  5.2  deg/sec 


would  be  expected.  The  yaw  rate  in  Figure  16  approaches  this,  but  is 
effected  by  limiting  in  the  telemetry  system. 

In  summary,  it  has  been  demonstrated  that  the  digital  implementation 
of  the  ZOH  controller  performed  properly  —  both  in  the  hybrid  simula¬ 
tion  and  in  VRA  flights. 
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a)  Theoretical  Response  of  Continuous  Controller 
to  Step  Input  of  Pedo!  (/  radian  Canard  Com¬ 
mand  Step) 


b)  /O  Hz  ZOH  Controller  with  Microprocessor 
Coupled  to  TR-  48  Analog  Computer  Response 
to  Redo!  Step  Input 


Figure  15.  Comparison  of  Theoretical  Continuous  Controllers  with  ZOH 

Controller  on  Hybrid  Simulator 


SECTION  IV 


SLEHER  CONTROLLER 


A.  INTRODUCTION 

This  section  describes  the  control  law  synthesis  procedure  and  vali¬ 
dation  process  for  the  slewer.  Because  of  "structural"  differences 
between  the  ZOH  and  slewer  first-order  equations  of  motion  for  the  YF-16 
it  was  necessary  to  "invent"  a  generalization  of  equivalent  stability 
derivative  (ESD)  model  matching  approach.  The  nature  of  the  solution  is 
described  first,  with  crucial  theoretical  proofs  relegated  to  appendi¬ 
ces.  A  low-order  example  is  given  to  illustrate  the  mathematical  pro¬ 
cess  and  a  general  four-step  synthesis  procedure  is  outlined. 

Next,  the  control  law  for  the  slewer/YF-16  combination  is  given, 
with  the  intermediate  numerical  details  given  in  an  appendix.  The  con¬ 
trol  law  is  theoretically  validated  by  running  step  responses  with  the 
theoretical  controller  working  against  an  analytical  model  of  the  open- 
loop  YF-16.  The  generalized  ESD  approach  does  an  excellent  job  of 
matching  the  simulator  to  the  model. 

Next,  the  details  of  validating  the  hardware  (microprocessor)  based 
control  law  against  an  analog  (TR-48)  model  of  the  YF-16  are  described. 
The  control  law  is  validated,  as  is  the  physical  implementation  of  the 
slewer  (in  terms  of  hardware  integrators  and  a  "software  differencing" 
format.)  Noise  problems  surfaced  in  the  next  stage  which  prevented  the 
flight  validation  phase  from  being  carried  out.  These  problems  are 
described  in  some  detail. 

B.  SYNTHESIS  USING  A  GENERALIZED  ESD  APPROACH 

A  radical  departure  in  the  format  of  the  first-order  form  occurs 
when  the  slewer  is  used  as  the  coupler.  The  first-order  form  (see 
Appendix  E  for  a  derivation)  becomes 
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^8k+l  m  ^8^8k  ^l^®k  +  ^^^ak-l  (69) 

while  model  equations  remain  the  same 

X%+1  "  $mxmic  +  ^m^m^  (?0) 

The  pseudo  inverse,  which  was  very  effective  in  the  ZOH  case,  will  no 
longer  work*  To  see  this,  transform  the  above  equations  and  substitute 
the  control  law 


5s  *  cF<Sm  +  cBxs  (71) 

Into  the  simulator  equation*  The  result  is  (see  Appendix  F): 

[iz  -  ♦g  -  (i^  +  r2z“1)cB]xg  -  [fj  +  z“1]cf6i 


Jm 


(72) 


[iz  “  $m]xm  “  rm6m 


Thus  an  exact  match  requires 


rc„i 

r  r 

F 

-  [rl  +  T2Z"1]"1 

m 

_  C]B  _ 

.  ♦m  ~  4>s  _ 

(73) 


which  Is,  In  general,  unattainable  for  two  reasons*  The  first  is  that 
[I*i  +  r2Z”l]  *  is  not  invertible*  The  second  problem  is  that,  even  if 
the  inverse  exists,  [Tj  +  ^z”!]  *  may  have  poles  exterior  to  the  unit 
circle  and  therefore  produce  an  unstable  controller.  Needed  is  a  syn¬ 
thesis  procedure  that  generalizes  the  ESD  approach,  can  treat  non-square 
matrices,  and  Insures  a  stable  controller.  The  Wiener-Hopf  approach 
satisfies  these  criteria  and  is  discussed  next* 
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C.  SOLUTION  VIA  VIENKR-HOPF 


Given  the  equation  set  (Equation  72)  the  W-H  equation  becomes  (see 
Appendix  F): 


(74) 


If  one  considers  [Cp  |  CB]  as  a  single  matrix,  and  thus  in  a  certain 
sense  one  unknown,  the  problem  posed  in  Equation  74  becomes  clear;  there 
is  only  one  (matrix)  equation  but  two  unknowns.  That  is,  treating 
[Cp  |  CB]  as  one  unknown,  the  other  unknown  is  the  matrix  ip.  All  we 
know  about  ip  is  that  it  must  have  poles  exterior  to  the  unit  circle; 
that  is,  a  property  of  ip  is  known.  The  trick  is  to  "pick"  a  [Cp  |  CB] 
matrix  that  is  stable,  substitute  it  in  Equation  74,  and  find,  after  the 
evaluation,  that  ip  has  the  desired  property. 


D.  SPECTRAL  FACTORIZATION  SOLUTION  METHOD 


Equation  74  has  the  form 


<pY  -  N  -  ip  (75) 

If  a  factorization  of  $  exists  such  that 

♦  *  $1*$1  (76) 

where  <p^  has  the  poles  interior  to  the  unit  circle  and  <p^  ■  <p'(l/z)  has 
poles  exterior  to  the  unit  circle,  then  the  solution  can  be  written  as 

Y  -  <h1[*liN]+  (77) 

where  £  ]+  signifies  a  partial  fraction  expansion  in  which  only  the 
terms  associated  with  poles  interior  to  the  unit  circle  are  retained. 

We  must  verify  that  Equation  77  produces  a  ip-vector  that  has  all  of  its 

poles  exterior  to  the  unit  circle: 


♦y  -  n  -  -  ♦ 

-  -  ♦uN]  -  * 


(78) 


Write 


♦Tin  -  [n3'N]+  +  [♦;1N]_  (79) 

where  [  ]+  Is  the  sum  of  the  partial  fraction  terms  interior  to  the  unit 
circle  and  [  ]  is  the  sum  of  the  terms  exterior  to  the  unit  circle. 
Thus 


♦i*u*iiN)+  -  (*;1n)+  -  u;1n)J  -  * 

and 


(80) 


(81) 


which  Is  guaranteed  to  have  all  poles  exterior  to  the  unit  circle, 
verifying  Equation  77  as  a  solution  of  the  W-H  equation. 

The  difficulty  with  the  approach  is  that  the  spectral  matrix  <t>  has 
to  be  factored  into  the  product  of  two  matrices,  one  which  is  analytic 
in  the  interior  of  the  unit  circle,  the  other  analytic  in  the  exterior 
of  the  unit  circle.  This  is  a  formidable  task.  However,  the  factoriza¬ 
tion  approach  does  give  key  insights  into  the  properties  of  our  particu¬ 
lar  W-H  equation,  a  matter  to  be  discussed  next. 

E.  KEY  FEATURES  OF  THE  SOLUTION 


From  the  previous  discussion,  if  the  spectral  matrix  can  be  written 


the 


♦  "  $  l 


det  <fr  *  det  <»>  det  $  i 


(82) 

(83) 
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and  we  see  that  the  poles  of  $  Interior  to  the  unit  circle  must  be  the 
poles  of  and  the  poles  of  $  exterior  to  the  unit  circle  must  be  the 
poles  of  ^  1  A  key  observation  is  therefore  that  the  poles  of  the 

solution  vector  Y  are  already  known;  we  need  only  determine  the  zeros 
of  Y.  To  do  this  assume  the  numerators  of  Y  are  unknown  polynomials 
divided  by  the  "stable"  poles  of  <j>.  Substituting  Y  into  Equation  75 
must  lead  to  a  ^  which  has  only  poles  exterior  to  the  unit  circle; 
therefore  the  numerator  polynomials  must  be  selected  so  that  the  poles 
of  <j>i  cancel  out  of  Equation  75.  This  is  the  "direct"  approach  for 
solving  W-H  equations  that  is  described  in  detail  in  Reference  3. 

F.  ILLUSTRATIVE  EXAMPLE 

Refocusing  our  attention  on  Equation  74,  we  show  how  the  thought 
process  described  earlier  can  be  used  to  effect  a  solution.  A  simple 
illustrative  example  is  used  to  facilitate  the  discussion.  Suppose 
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rl  -  r2*  ” 


z  +  1  .75z 


-• 25z  +  1 


(88) 


The  spectral  matrix  *  Is 


*  -  [rI  +  r2zJ[ri  +  r^-1] 


[ 


z2+2. 5625z+l 
•  8125z+l. 75 


25z 


1.75z2+.8125z 

2+2.0625z~.25z-.25j 


(89) 


det  *  -  -.25z4  +  1.0625z2  -  .25  -  -25[z4  -  4.25z2  +  l] 


(90) 


The  roots  of  Equation  90  are 


D  -  (z  -  .5)  (z  +  .5)  (z  -  2)  (z  +  2) 


and  therefore  select 


(z  -  • 5) (z  +  .5)  -  z*  -  .25 


(91  ) 


as  the  poles  for  all  the  entries  of  Cp  end  Cg .  Equation  74  has  the 
form,  for  Cps 


[ 


*2-2.5625z+1 
. 8125z+l. 75 


1. 75z2+.  8125z  J  fao*29e1z+a2  bQZ2-*^  jz+b2 
-.25z2+2.0625z-.25J  Lc( 


co*2+ciz-H:2  dQZ2+diz+d2  J 


(z  -  .5)(z  +  .5) 


[ 


2. 5z+l  2.75z+2" 
-• 5z+3  -. 25z+3. 


z*(z  -  .5)  (z  +  .5) 


z  (z  -  .  5)  (z  +  .  5)  - 


(92) 
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In  Equation  92  the  (1*1  has  been  augmented  by  a  zA  factor,  in 

order  to  force  a  common  denominator.  We  must  pick  the  unknown  coeffi- 
dents  of  Cp  so  that  each  numerator  of  Equation  92  contains  z  (z  -  .5) 
(z  +  .5),  in  order  that  ^  will  be  free  of  poles  interior  to  the  unit 
circle*  One  property  of  the  solution  is  therefore  immediately  apparent; 
setting  z  «■  0  gives,  for  the  numerator  of  Equation  92, 


1 

o 

•H 

1 _ 

a2  b2 

1.75  .25 

C2  d2 

(93) 


Since  det 


1 

_1. 75 


we  are 


assured  there  is  the  unique  solution 


a2 

c2 


(94) 


This  is  a  general  result;  all  numerator  entries  of  Cp  and  CB  must  con¬ 
tain  a  free  z.  The  free  z  in  Cp  and  Cg  will  force  a  cancellation  of  one 
of  the  z's  in  the  first  term  of  Equation  92.  This  step  is  of  such 
Importance  that  we  emphasize  it  by  simplifying  Equation  92: 


z2+2.5625z+l 

1.75z2+.8125z 

agz+ai 

bgz+bi 

. 8125z+l. 75 

-. 25z2+2. 0625z-.  25 

_c0z+ci 

dQz+d 1. 

z  (z-.5)(z+.5) 


2. 5z+l  2.75z+2 


. 5z+3  . 25z+3  J 

-  (z  -  . 5) (z  +  .5)  -  *  (95) 

z (z-. 5) (z+.5) 
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There  is  now  only  one  free  z  in  the  denominator  of  Equation  95.  The 
numerators  of  Equation  95  must  all  contain  a  free  z,  in  order  for  the 
z  in  the  denominator  to  cancel  out  of  tju  Setting  z  -  0  gives 


1  0 

1.75  -.25 


1 

*— 1 
XI 

<0 

1 _ 

"l  2' 

- 

_C1  dl. 

3  3 

(—  5)  ( .5) 


(96) 


or 


H 

f0 

1 _ 

bl' 

'  —  25 

•e 

—  5 

m 

„C1 

dl. 

1.25 

-.5 

(97) 


The  next  step  is  to  pick  (Sq,  bg,  Cq,  dp)  so  that  A  -  (z  -  .5)(z  +  .5) 
is  not  a  root  of  <j*.  However,  for  those  values  of  z  such  that  A(z)  -  0, 
Equation  95  simplified  tremendously  and  we  obtain  another  general  result 


$Cp  s  0  when  A(z)  -  0 


(98) 


This  result,  however,  must  be  treated  with  care;  it  does  not  mean  that 
Cp  ■  0  for  these  values  of  z  since  does  not  exist.  That  is. 


det  =  0  for  z  such  that  A(z)  -  0 

The  correct  usage  of  Equation  98  is 


(99) 


[♦] 


8qZ  +  a  L  bQz  +  bx 


cqz  +  cj  dQZ  +  dj 


-  0 


(100) 


’a0z  b0z" 

"  8 1  bl" 

[♦] 

m  — <p 

CQZ  dQZ 

_  m 

_C1  dl. 

(101) 
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Clearly,  since  det  ♦  -  0  for  z  -  0.5,  -0.5,  the  two  equations  of  $  are 
linearly  dependent;  thus  we  need  only  one  of  them.  Suppose  we  choose 
the  first  row: 


[(z2  +  2.5625z  +  1)  |  1. 75z2  +  .8125*] 


aQz  +  aj  bQz  +  bj 


Cqz  +  cj  dgz  +  dj 
where  z  ■>  .5,  -.5 


When  z  -  0.5, 


(102) 


‘.5a0 

+  a^  .  5bg  +  b 

.84375] 

_.5c0 

+  c^  •  5dg  +  d 

(a0»  bo* 

Cq,  dg)  in  terms  of 

r-z2 

-  •  25z  z2  - 

_ 

L  2z2 

+  1.25*  z2  - 

(z  -  .5) (z  +  .5) 


(103) 


(104) 


Using  the  same  procedure,  except  substituting  <|>m  -  $s  for  Tm>  gives 


CB 


-. 25z  z2  -  .5z 

. z2  +  • 25z  z2  -  . 5z 
(z  -  . 5) (z  +  .5) 


(105) 


This  Illustrative  example  was  chosen  so  that  +  I^z-*  would  be  invert¬ 
ible  and  give  the  reader  an  opportunity  to  verify  the  answer  using 
Equation  83.  However,  bear  in  mind  that  the  W-H  approach  works  when 
1*1  +  r^z'^ia  not  Invertible. 


'  -- 


Before,  leaving  the  example,  perform  the  computation 


cplcB  “ 


1/3  1 


(106) 


and  observe  another  property  of  the  solution,  namely  that  Cp  Cg  is 
always  a  matrix  of  constants.  It  is  given  explicitly  by  the  equation 
(see  Appendix  G  for  proof): 


Cf^b  -  (rlrm)“lrl(^m  “  $s) 


(107) 


G.  GMTKRAI.  PROCEDURE 

The  steps  used  in  the  illustrative  example  are  a  good  basis  for  a 
general  procedure.  First,  let 


Cp  * 


(108) 


where  Cp^,  Cg*^  are  the  polynomial  adjoint  matrices  of  A,  and  rewrite 
Equation  74  as 


[ri  +  r2z][riz  +  r2)[cfdJ  ;  cgdJ]  -  [T1  +  r2z][rm(^m  -  <frs)]A 


There  are  four  steps  in  the  procedure: 

1)  Evaluate  det[(Ti  +  T2z)(riz  +  ) ] ,  select  the 

interior  roots  for  A(s). 

2)  Set  up  Cp^  and  Cg*^  with  unknown  polynomial 
entries.  The  order  of  each  polynomial  is  equal 
to  the  order  of  A(z),  with  all  z  terms  set  to 


(109) 


[_  *  w  m  . 


3)  Solve  for  Cpd^(0),  Og^CO)  using 

[Cfdj  |  cSdJ]  -  (rlr2)'1  ri[rn](^  -  ^)a(O)  (lio) 

4)  Find  Che  remaining  equations  by  evaluating  appro¬ 
priate  rows  of 

(rI  +  r2z)(riz  +  r2)[cfdj  j  c|dJ]  -  o  ,  a(z)  -  o  (ill) 

Next,  we  apply  the  procedure  to  the  slewer  controller  ise  and  tabulate 

the  Cp  and  Cg  matrices. 

H.  SLEWER  CONTROL  LAW 

Figure  17a  has  a  matrix  Cp  premultiplying  the  Cp  matrix.  As  in  the 
case  of  the  ZOH,  this  is  the  matrix  which  couples  the  model  inputs  to 
force  a  wings-level  turn.  As  in  the  case  of  the  ZOH,  it  is  a  3  x  1 
matrix,  but  the  entries  are  slightly  different. 


Figure  17.  Block  Diagram  for  Software  Realization 
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The  idealized  configuration  of  Figure  17a  must  be  modified  for  the 
alever  alnce  we  have  no  phyalcal  implementation  of  a  completely  hardware 
nature  for  implementing 


M 


fl  -  e-sTl2 

Ta2 


(112) 


We  choo8e  to  Implement  the  alever  by  "breaking  up"  Equation  112  accord¬ 
ing  to  Figure  17b.  In  Figure  17b  the  basic  transfer  function  of  the 
dlgltal/analog  converter  as  a  ZOH  has  been  retained,  while  the  addi- 
tional  (1  -  e“B  )/T  has  been  moved  back  into  software  as  (z  -  1)/Tz.  A 
physical  (analog)  Integrator  has  been  inserted  as  a  prefilter  to  the 
controllers.  This  is  an  implementation  that  theoretically  matches  the 
characteristics  of  Equation  112.  The  theoretical  details  are  discussed 
in  Appendix  H. 

We  will  solve  for  the  CF  and  Cg  matrices  using  the  procedure  out¬ 
lined  above.  The  conqmtatlonal  details  are  very  extensive  and  for  that 
reason  we  relegate  the  input  data  and  numerical  computations  to  Appendix 
I.  However,  we  must  first  discuss  the  structure  of  the  control  law. 

Each  entry  of  CB  and  CF  is  a  ratio  of  polynomials  in  z.  Dividing 

through  by  the  highest  power  of  z  gives  a  large  number  of  terms  that 

-I  -2 

will  be  multiplied  by  delay  operators  such  as  z  ,  z  ,  etc.  We  will 
group  all  such  operations  involving  "past"  values  of  the  data  in  a  back¬ 
ground  computation.  We  proceed  as  follows,  defining: 


z[B2Z^  +  Bjz  +  Bg] 


A(z) 


CpCp  ” 


z[C2z^  +  Cjz  +  Cq] 


A(z) 


(113) 


(114) 


From  Figure  17b, 


6  ■  Cp6p  +  CjjX^ 


z [c2z^  +  Cq]  j  z[b2z^  +  B^z  +  Bq] 

- Iw - Sp  +  " 


A(z) 


Dividing  through  by  zJ  gives 


[l  +  a2z“l  +  ajz“^  +  aQZ“3Jfi  ■  (c2  +  Cjz- ^  +  Cgz~2)5T 


+  (b2  +  B12-1  +  B0z"2)xT  (117) 


The  recursion  equation  is 


6k  “  ~a26k-l  “  *1^-2  “  a06k-3  +  c26Pk  +  cl6Pk-1  +  c06Pk_2 


Define  the  "background”  computation  as: 


bk  ”  -«z5k-l  -  al^k-2  -  a0fik-3  +  Cl6Pk_i  +  C06Pk_2 


+  Blxk-1  +  ®0xk-2 


Thus  the  "foreground"  Is 


6k  -  C25p,  +  B2Xk  +  Bk 


In  addition,  we  need 


{2  -  i  t«k  -  Sk-J 
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;  ;v  v; .  v , :  v 


(116) 


+  B2Xk  +  BiXk_i  +  BgXk_2  (118) 


(119) 


(120) 


(121) 


when  new  data  are  taken  In,  the  microprocessor  will  compute  the 
first  component  of  Equations  118  and  121  (the  rudder  control  law), 
enable  an  interrupt  and  output  to  the  rudder.  It  will  then  repeat  this 
process  for  the  side  force  controller  and  the  aileron  controller.  In 
the  remainder  of  the  frame  time,  the  background  computation.  Equa¬ 
tion  119,  will  be  updated.  The  data  for  the  matrices  are  given  in 
Table  5. 


I.  VALIDATION:  THEORETICAL  CONTROLLER,  THEORETICAL  YF-16 

The  control  law  defined  by  Equations  118,  119,  and  120  was  theoreti¬ 
cally  validated  by  running  step  responses.  These  are  shown  in  Fig¬ 
ure  18.  The  agreement  between  the  slewer-based  system  and  the  contin¬ 
uous  baseline  is  very  good.  The  step  responses  are  also  tabulated  in 
Table  1-1  of  Appendix  H  (these  can  be  compared  against  Table  4  of  Sec¬ 
tion  III). 

J.  VALIDATION:  HARDWARE-EASED  SLEWER,  ANALOG  YF-16 

Figure  19  shows  a  time  response  from  a  hybrid  simulation  of  the 
alewer-controlled  YF-16.  The  control  law  is  implemented  in  the  micro¬ 
processor  with  the  YF-16  model  on  the  TR-48  analog  computer. 

The  input  to  the  YF-16  model  is  a  step  input  of  direct  force  command 
(pedal)  from  a  switch  on  the  TR-48.  This  analog  signal  enters  the 
microprocessor  and  Is  processed  according  to  the  slewer  control  law, 
which  puts  out  canard,  rudder,  and  aileron  surface  rate  signals. 

The  YF-16  canard  rate  command,  is  shown  in  the  top  trace  of  Fig¬ 
ure  19.  There  is  a  "ringing"  noise  signal  superimposed  on  the  rate  com¬ 
mand  signal.  This  noise  is  much  larger  than  the  rate  signal  itself. 
This  suggests  an  improper  scaling  of  one  of  the  channels  which  is  allow¬ 
ing  a  noise  response  of  an  element  of  the  matrix  digital  filter  to  be 
Incorrectly  added  to  the  pedal  command.  We  shall  have  more  to  say  about 
this  noise  problem  at  a  later  point. 
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The  Cp  matrix  of  Figure  17  breaks  the  signal  up  into  three  compo¬ 
nents  to  generate  three  surface  rate  signals  which  are  Integrated  by  an 
analog  integrator  on  the  TR-48.  This  simulates  the  analog  integrators 
in  the  VRA  "mixing  box."  The  YF-16  canard  position  command,  from  its 

e 

integrator,  is  shown  in  the  second  trace.  It  is  seen  that  the  6q  signal 
has  been  properly  Integrated  to  a  step,  with  a  small  10  Hz  noise  compo¬ 
nent  superimposed.  The  slow  drift  in  <5C  (and  the  response  signals)  is 
due  to  the  bias  in  the  D/A  converters,  which  input  the  analog  integra¬ 
tors.  This  bias  could  not  be  completely  eliminated  by  the  integritor 
balancing  circuits. 

The  yaw  rate  response  shown  in  the  bottom  trace  has  a  "steady-state" 
value  of 


0.08  (deg/sec) /deg 


which  compares  well  with  the  calculated  value  of  0.059  (deg/sec) /deg 
(considering  the  drift).  The  transient  response  shows  the  appropriate 
rise  time  and  overshoot. 

The  sideslip,  roll  rate,  and  bank  angle  traces  remain  essentially 
zero  (when  the  integrator  drift  is  allowed  for)  as  required  for  the 
wings-level  coordinated  turn. 

It  was  necessary  to  keep  signal  levels  low  because  larger  Inputs 
caused  the  noise  component  to  limit  the  surface  rate  channels.  However, 
because  the  simulation  is  linear,  it  can  be  seen  that  the  slewer  con¬ 
troller  performs  properly  over  the  frequency  region  of  interest  for  con¬ 
trol  of  the  YF-16  dynamics. 

In  essence,  the  problems  noted  in  the  traces,  such  as  the  unusually 
sensitive  response  to  digital  noise,  are  not  properties  of  the  slewer 
control  law  but  properties  of  any  recursive  control  law.  That  is,  D/A 
and  A/D  interfaces  imply  a  ratloing  between  physical  variables  and 
machine  variables.  If  the  gain  of  a  machine  variable  is  too  high,  the 
system  will  always  be  responding  to  digital  noise  or  bit  dither.  Thus 
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the  problem  is  to  properly  proportion  physical  variables  and  machine 
variables,  a  process  which  requires  a  systematic  technique  for  perform¬ 
ing  the  scaling.  Unfortunately,  our  original  plan  for  gathering  the 
data  needed  for  the  scaling  was  not  workable.  The  microprocessor  speed 

prevented  a  real-time  printout  of  the  states  and  controllers,  since  we 
were  working  with  a  TR-48,  which  would  only  work  in  real  time.  These 
problem  are  described  next. 

K.  GAIN  ALLOCATION  AND  NOISE  PROBLEMS 

One  of  the  major  obstacles  to  successful  evaluation  of  the  slewer 
control  law  and,  in  fact,  a  major  consideration  in  any  digital  control 
environment,  may  be  termed  the  gain  allocation  problem.  Consider,  for 
Illustration,  the  zero-order  hold  case  (Figure  20a) . 


3x3 


a) 


b) 


Figure  20.  Zero-Order  Hold  Gain  Configurations,  a)  Standard 
Nomenclature,  b)  Equivalent  Configuration  vs.  Implemented 
in  DFCS  Control  Hardware/Software 


In  the  present  digital  implementation  the  computations  are  struc¬ 
tured  according  to  the  alternative  form  (Figure  20b). 

If,  for  checkout  purposes  and  computational  convenience,  one  desires 
to  enter  the  gains  for  Cp  and  Cg  in  physical  units,  it  is  necessary  to 
introduce  additional  scale  vectors  (Figure  21),  where  A  and  B  are  vec¬ 
tors  of  dimension  3  and  4,  respectively,  which  convert  command  and  state 
A/D  voltages  to  physical  units;  and  C  is  a  vector  of  dimension  3,  which 
converts  control  and  physical  units  to  the  correct  voltage  units  for 
driving  the  servo  buffer  amplifiers. 

This  particular  form  worked  well  in  earlier  assembly  code  versions 
of  the  control  program.  However,  in  the  later,  more  computationally 
advanced  versions,  performing  the  vector-matrix  operations 

6X  x  B  *  Cg  x  C  and  6m  x  A  x  Cp  x  C 

could  not  be  reliably  accomplished  within  the  100  msec  time  frame.  This 

problem  was  solved  by  premultiplying  the  scaling  vectors  prior  to  enter¬ 
ing  the  control  loop,  thus: 

A  A 

Cp  -  A  x  Cp  x  C  and  Cg  =  B  x  Cg  x  c 

Now  the  real-time  calculations  revert  to  the  simplified  form  (Figure  22) 
and  are  easily  accomplished  within  the  designated  time  frame.  Although 

the  slewer  case  is  intrinsically  more  complex,  the  same  operational 

principles  apply  to  the  digital  code. 

A  A 

In  the  ZOH  case  th^  magnitudes  of  Cp  and  Cg  fell  between  ±1.2.  In 
the  slewer  case,  however,  the  presence  of  the  hardware  integrators 

AAA 

forced  the  conditioned  gains  (Cg,  Cj ,  C2)  to  vary  over  a  much  wider 
range  (±40).  While  the  range  of  intermediate  products  was  not  observed, 
we  estimate  that  transient  values  exceeding  10  may  have  occurred  during 
slewer  control  law  operations.  The  presence  of  high  gains  and  uneven 
gain  distributions  makes  digital  systems  highly  vulnerable  to  digital 
artifacts  such  as  bit  dither,  roundoff,  and  truncation  errors.  In  the 
slewer  case,  additional  problems  of  analog  noise  sensitivity  in  the  A/D 
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Figure  21.  DFCS  Gain  Configurations  Illustrating  Within-Loop 
Placement  of  Scale  Vectors 


3x3 


Figure  22.  Resultant  Gain  Configuration  Illustrating  Simplification 
Brought  About  by  Pre-Multiplication  of  Matrices 
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converters  and  offset  sensitivity  in  the  hardware  integrators  compounded 
our  difficulties* 

It  should  be  noted  here  that  sensitivity  problems  of  this  nature  are 
extremely  hard  to  identify,  especially  due  to  the  dynamic  requirements 
of  the  slower*  It  is  nearly  Impossible  to  capture  intermediate  values 
with  the  slower  running  in  real  time  due  to  the  intolerance  of  the  frame 
epoch  for  additional  code,  especially  I/O.  Another  approach  would  be  to 
slow  down  the  TR-48  to  match  the  frame  rate  of  a  modified  slewer  with 
I/O  in  the  control  loop*  This  procedure,  however,  is  resisted  by  opera¬ 
tional  limitations  of  the  TR-48.  A  general  solution  of  minimizing  the 

AAA 

range  of  conditioned  gains  (Cg,  ,  C2)  was  attempted  by  using  Internal 
analog  scaling  to  reduce  the  internal  gain  distribution  from  +40  to  +9* 

This  was  somewhat  successful  in  that  it  resulted  in  the  only  stable  data 
obtained  on  the  slewer.  However,  analog  noise  external  to  the  digital 
system  was  exacerbated. 

Future  efforts  would  definitely  include  a  systematic  attack  on  scal¬ 
ing  problems  of  this  type,  since  they  appear  to  be  the  root  of  many 
digital  control  system  problems. 

L.  THE  GRACEFUL  ENTRY  PROBLEM 

Although  it  is  generally  considered  appropriate  to  trim  an  aircraft 
to  stable  level  flight  prior  to  engaging  an  autopilot,  it  has  long  been 
recognized  that  some  sort  of  autotrlm  mechanism  is  necessary  for  smooth 
engagement  since  perfect  level  flight  is  seldom  extant  at  engage  time. 

The  Princeton  VRA  represents  no  exception  to  this  rule.  Minimiza¬ 
tion  of  engage  transients  is  extremely  critical  to  the  Navlon  due  to  the 
high  angular  accelerations  available  at  the  control  servos. 

The  VRA's  analog  system  has  an  active  continuous  autotrim  which  is 
highly  effective.  For  the  micro-DFCS,  however,  it  is  necessary  to 
Implement  the  autotrim  in  software.  The  Figure  23  flow  diagram  illus¬ 
trates  how  this  is  done,  again,  using  the  zero-order  hold  case  as  an 
example.  An  important  consideration  in  dealing  with  control  law  calcu¬ 
lations  on  the  micro-DFCS  is  that  the  numerical  basis  for  the  A/D  and 
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Figure  23.  ZOH  Auto-Trim  Algorithm 


D/A  peripherals  Is  offset  binary,  l.e.,  0  -  -10V,  2048^  -  211  -  OV, 
4096jg  ■  2*^  ■  +10V.  Since  the  same  numerical  basis  is  used  for  both 
A/D  and  D/A  operations  It  Is  considerably  more  efficient  to  perform 
control  law  calculations  with  the  offset  in  place.  This  means,  using 
the  20H  case  as  an  example,  for  small  state  and  command  perturbations 
centering  around  0  db  2.0  VDC,  control  law  products  are  in  the  range 
1.0  x  10^  <  P  <  5.0  x  10^  or  11-13  bits.  Since  the  autotrim  calcula¬ 
tions  which  are  applied  to  the  Z0H  case  operate  in  the  same  magnitude 
range  they  pose  no  threat  of  digital  saturation  even  with  the  binary 
offset  present. 

In  the  slewer  case,  however,  an  entirely  different  set  of  parameters 
was  found  to  apply.  First,  there  is  an  order  of  magnitude  differential 
between  A/D  incoming  volts  and  D/A  output  volts  due  to  the  present  of 
the  hardware  integrators.  Second,  the  presence  of  the  11  bit  binary 
offset,  coupled  with  the  relatively  larger  gains  and  greater  number  of 
computational  terms  produced  intermediate  products  which  apparently 
caused  saturation  of  the  9511  match  chip.  For  these  reasons,  the  ZOH 
autotrlm  algorithm  would  not  work  for  the  slewer.  Additionally,  it  was 
found  that  the  slewer  control  loop  computations  were  numerically  incor¬ 
rect  when  performed  in  offset  binary.  Although  modifying  the  control 
loop  routine  to  take  out  the  offset  immediately  following  A/D  input,  and 
reinserting  it  immediately  prior  to  D/A  turned  out  to  be  fairly  simple, 
modifying  the  autotrim  required  a  major  code  rework  which  could  not  be 
accomplished  in  the  remaining  time  at  Princeton. 

The  slewer  was  evaluated  on  the  TR-48  by  initializing  the  analog 
computer  in  a  manner  that  simulated  autotrim.  (Digital  autotrim  was 
disabled  for  these  test  runs.)  An  analogous  simulation  in  the  VRA  could 
not  be  accomplished  due  to  the  design  of  the  VRA  analog  system.  This 
was  the  major  impediment  to  flight  testing  the  slewer  on  the  VRA. 


SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

From  a  theoretical  viewpoint,  both  the  direct  synthesis  ZOH  and 
8lewer  control  laws  achieve  the  design  objective  of  emulating  the  base¬ 
line  wings-level-turn  direct  side  force  mode.  In  the  ZOH  case  the 
equivalent  stability  derivative  (ESD)  model  follower  simulator  approach 
was  enhanced  by  the  innovative  use  of  a  pseudo  Inverse. 

A  major  theoretical  addition  to  in-flight  simulator  technology  was 
achieved  in  the  slewer  design  with  the  introduction  of  a  Wiener-Hopf 
optimization  algorithm  for  augmenting  the  basic  ESD  approach.  This  is 
of  great  value  in  that  other  configurations  (other  than  the  slewer)  lead 
to  the  same  class  of  problem.  For  example,  digital  controllers  utiliz¬ 
ing  ZOHs,  wherein  computer  throughput  delay  is  of  concern,  produce  the 
exact  same  problem  formulation. 

The  ZOH  control  law  was  hardware  validated  against  both  the  TR-48 
(functioning  as  a  surrogate  YF-16)  and  the  Navion  VRA  (with  analog 
closures  to  simulate  an  open-loop  YF-16). 

The  slewer  controller  was  validated  in  hardware  against  a  surrogate 
YF-16  (TR-48).  This  is  significant  in  that  the  input  to  the  actuators, 
from  the  hardware  Integrators,  possessed  all  the  correct  attributes  even 
though  a  "continuous"  slewer  term,  1  -  e”sT,  was  implemented  in  software 
as  (z  -  1)/Tz.  Thus,  one  possible  physical  realization  of  the  slewer 
coupler  was  realized. 

On  the  negative  side,  largely  unanticipated  problems  involving  the 
scaling  of  physical  and  machine  variables  between  the  A/D  computer  and 
D/A  modules  led  to  an  inordinate  and  unacceptable  level  of  digital  noise 
(bit  dither,  rounding,  etc.).  This  noise  was  a  major  impediment  to 
flight  validation  of  the  slewer.  Another  major  obstacle  was  the  lack  of 
an  adequate  digital  autotrlm  algorithm  for  the  slewer.  This  created  an 
engage-transient  or  "graceful  entry"  problem  which  precluded  flight  test 
of  the  slewer. 

B.  KECOMfENDATIONS 

There  are  four  recommendations.  First,  two  software  routines  need 
to  be  written  for  the  slewer/microprocesror /Navion  combination: 
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1)  An  autotrim  algorithm  for  the  YF-16  program  that 
provides  graceful  entry  into  the  test  run. 

2)  A  program  for  optimizing  the  relative  ratios 
between  physical  variables  and  machine  variables 
so  as  to  eliminate  or  minimize  the  digital  noise 
problem. 

With  the  control  of  the  digital  noise,  the  integrator  drift  problem  will 
probably  disappear.  If  it  remains,  then  we  need: 

3)  A  slewer  realization  that  resets  the  integrator 
every  T  seconds. 

Implementation  of  Recommendations  1  and  2  should  be  straightforward 
since  both  Issues  underwent  preliminary  attack  immediately  prior  to  the 
termination  of  our  last  flight  test  effort  at  Princeton.  The  autotrim 
problem  is  easily  solved,  but  was  unanticipated  until  our  last  minute 
review  of  slewer  computations  was  made.  The  realization  of  the  first 
two  recommendations  (and  if  necessary  the  third)  would  set  the  stage  for 
carrying  out  the  fourth  recommendation: 

4)  Pilot  evaluation  flights  to  furnish  paired  com¬ 
parisons  between  the  ZOH,  slewer,  and  baseline 
analog  configurations. 
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APPENDIX  A 


DAXA  FOR  FIRST-ORDER  FORMS 


This  appendix  gives  data  for  the  first-order  forms  of  the  equations 
of  motion .  The  continuous  system  has  the  form 

X  -  (AMAT)X  +  (BMAT ) 6 

whereas  the  ZOH  forms  are 

%+1  -  (FMAT)Xk  +  (GMAT)6k 

Data  are  for  open-loop  Navi  on,  YF-16  CCV  with  canard  (the  bare  air¬ 
frame),  and  the  YF-16  CCV  with  canard  and  FCS  loop  in  (the  closed-loop, 
reduced-order  model).  Sample  rate  varies  from  1  to  100  Hz. 
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CASE  II  —  T«1 


YF— 16  CCV  Canard  FCS  Loops  Closed 
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_. 208547660 8E -02 
.2256656150 
-.  1272278461 
.3927491796 
2.201845451 


.37950000006-01 

-25.00000000 

e. 


-.2yj!5M7y:?p-ei 
-.26549196935-01 
-.53072417236-01 
_ .5551S06333E-01 


72 


NAVION 


CASE  I  —  T-.2 


AHAT(ROU 

1)» 

-.7779000000 

4.680000000 

-.4320000000 

0. 

AHAT  ( ROU 

2)* 

“1.000060000  - 

-.3S56000000 

0. 

. 1720000000 

AHAT ( ROU 

3  >• 

1.270000000 

-12.80090000 

-6.600090000 

0. 

AHAT ( ROU 

4  )• 

It  ______ _ _ 

0. 

1.000000000 

0. 

tiwmfotr 

1  >• 

-6. 109099900 

~  2.419000009 

=731X00 0 0000 - 

BHAKROU 

2  )• 

•72S0000000E-01 

.2370000000 

e. 

BHAT(ROU 

3)- 

.7700000000 

0. 

21.00000000 

_bhat<rou 

9. 

0. 

0. 

iJ- 

.7638201656 

.8734526537 

-. 4088771668E-01 

. 15S4957149E-01 

FUAT (ROU 

2  )• 

-.1720708623 

.8410590820 

• 7465474683E-02 

. 32 15391975E-01 

FrATlROU 

3)* 

.2775824494 

-1.226753688 

.2575694151 

-.2739705798E-01 

FHAT(ROU 

4)- 

•2763312714E-01 

-.1592852208 

.1103686365 

.9979399271 

“CMftTTRO'tr 

1>- 

-1.088635063 

.4524869698 

-.1685774105 

CHAT (ROU 

2)* 

.1251228024 

.  . 3970933440E-O3 

. 17S3197902E-01 

CHAT (ROU 

3  )• 

-. 95 1264039 IE-01 

.2884523906E-01 

2.369064565 

CHAT (ROU 

4  )• 

-.1 6507852 17E-02 

. 1 560 46756 IE-02 

.2817568625 

YF-16  CCV  Canard 

Bare  A/F 

CASE  III  - 

-  T» 

.2 

AHAT (ROU 

I)- 

-.4310900009 

10.20000000 

-.4160000000E-01 

0. 

Af1AT(ROU 

2)- 

-1.000000009 

-.3060000000 

0. 

. 3880000000E-0 1 

AHAT (ROU 

3  )• 

1.670000000 

-SO. 00000000 

-2.330000000 

0. 

AHAT (ROU 

4)- 

0. 

9. 

1.000000000 

0. 

inAYTKOU 

TT?- 

“47T7W0BMS - 

- 5C9295W000 

27170000900 

BOAT (ROU 

2)» 

.318e00OO00E-01 

.19100009O0E-O1 

-.3570000000E-O1 

BHAT(ROU 

3  )• 

7.630000009 

5. 930000090 

49.10000000 

BHAT(ROU  _ 

4  )• 

0. 

0 

A. 

FHAT(ROU 

1  )• 

.7318284667 

1.798616821 

-.542834G869E-02 

. 73720 47226E-02 

FHATCROU 

2)- 

-.1730440775 

.7523851752 

.  1287590071E-02 

.  70276367 13E-0H 

FHAT(ROU 

3  )• 

1.024821400 

-6.888606229 

.6225241588 

-.3081227258E-01 

FHATCROU 

4  )• 

.8309475710E-<H_ 

-.7941307365 

.1595835351 

.9978105376 

^SflArfROU- 

nr 

- . 7447458733 

' .8261207402 

.3495396442 

CHAT (ROU 

2)« 

.8316230277E-01 

.-.8098278965E-01 

- . 4 190035662E-O1 

CHAT (ROU 

3  )• 

.8458638785 

1.315060244 

8.044217665 

CHAT (ROU 

4  >• 

. 1094266806 

.  1232560490 

.8575517720 

YF- 

16  CCV  Canard  FCS 

Loops  Closed 

CASE  II  - 

-  T- 

.2 

- 

AHAT(RCU 

1  )• 

-3.105000000 

8.917000000 

-.2720000000 

0. 

AHAT (ROU 

2  )• 

-.9796000000 

-.2965000000 

0 . 

. 3795000009E-01 

AHAT (ROU 

3  )• 

6.564000000 

-47.71000000 

-8.359000000 

-25.00000000 

AHAT (ROU 

4). 

0. 

0. 

1.000000000 

0. 

■"TTTATTROir 

1  )  • 

—=37925000000 

- 47765000000 

1.897000000 

BHATI ROU 

2>- 

.298P00PO00E-O1 

.  18O2eO0OO0E-0i 

- . 336O000000E-O1 

BHAT(ROU 

3)« 

7.178000000 

5.655000000 

49 . 60O00C00 

BHAT ( ROU 

4  >• 

0. 

0 .  ,  „  , 

0  • 

“FttfiTTROUT 

"TPT 

— rwrsrsm - 

- r730902S3T4 

-.1275477179E-01 

•6458359475E-01 

FHAT(ROU 

2)* 

-.1314973514 

.7988423546 

. 26748 45S25E-02 

. 1 37373 1217E-02 

FHAT(ROU 

3)« 

.7113296577 

-2.922093678 

. 1 131844369E-01 

-2.047166570 

FUAT (ROU 

4»- 

•9096334565E-01 

-.4705547884 

. 8 1 17236063E-01 

.7145802362 

tJfTATTROlT" 

Tpr 

-.5585048934 

652327 1267 

71419336178 

CHAT (ROU 

2  >• 

.66921 428 42E-01 

-.6804147674E-01 

-. 2395987392E-01 

CHAT (ROU 

3  )• 

.2115638959 

.8840816074 

4.214517195 

CHAT (ROU 

4J« 

.5394908806E-01 

•9566196556E-01 

.5774663642 

73 


NAVION 


CASE-T-.l 


amatirou  i >•  -.7770000000 

AflAT(ROU  2>«  -1.000000000 

>r?AT(ROU  3)-  1.270000000 

lltATtROU  4>«  e. 

ffwsTTWOO  i >  •  ~=t ;  rwwgrero - 

BnAT(ROU  2)-  .725000O000E-01 

BrtAT(ROU  3)*  .7700000000 

EMATCROU  4>«  0. 

“TWirrrROi/TTi — 773053752  si 

F»1AT(R0U  2)-  -.9363477977E-01 

FrtAT(ROU  3 )•  .1378664855 

,  -FHAT(ROU_  4>»  .6760554231E-02 

GflAT(ROU  1 )•  -.5814275833 

GrtAT(ROU  25-  .3637102593E-01 

GflAT  ( ROU  31«  .  1 141 512336E-01 

GflAT  (ROU  4).  . 1636999808E-02 


4.680000000 

-.3556000000 

-12.80000000 

0.  _ 

~ 2!74T0030000 
.2370000000 
0. 

0.  _ 

74598273385 
.9418349233 
-.8863597096 
-.5026331 29 4E-01 
.2353585281 
.  1 154709769E-01 
.4380530530E-02 
.  1 147306212E-03 


-.4320000000 

0. 

-6.600000000 
1.000000000 
"=73  740000000 
0. 

21.00000000 

0.  _  _ 
“2987563 1 93E^T 
. 2364 1 84296E-02 
.5142138486 
.731 150S762E^01_ 
- .6556590545E-01 
. 3265167507E-02 
1.533293396 
.8514173945E-01 


.  1720000000 


• 3990407457E-02 
. 167G410940E-01 
- .8645289826E-02 
.9996937882 _ 


CASE  111  —  T-.l 


YF-16  CCV  Canard  Bare  A/F 


ArtAT (ROU  1 »•  -.4310000000 

Art at ( rou  2 )•  -l.eeeeooeea 

ArtAT  <  ROU  3)-  1.670000000 

ArtATCROW  4>«  0. _ 

BrtATTROU  17=  -4 . 1V000O000 

SHAT (ROU  2)>  .318O000O0OE-01 

BflATCROU  31-  7.630000000 

BflAT  <  ROU  4  )■  0. _ 

"VTTOTOfOU  1  )•  7908561 1123 

FflAT  (ROU  2»-  -.9471662993E-01 

FrtATCROU  3)-  .3670070570 

FflAT ( ROU  4  >•  .  152 4291 762E-01 

TSWT  (ROU- lT-“'-‘.40 1 A594654 - 


WAT  (ROU  ~ IT.-' -‘.40 1 0594654 - 

CHAT (ROU  2)-  .2334360234E-01 

GflAT  ( ROU  3)-  .6097314035 

OAATCROU  41-  .3320928S55E-01 


10.20000000 

-.3060000000 

-50.00000000 

0. 

476200 00000 - 

. 191 O00000OE-0! 
5.930000000 
0. 

79755276552” - 

.9204149530 

-4.234745790 

-.2246473026 

“4443513628 - 

-. 20423525 3 IE-01 
.S949648390 
. 2942249255E-01 


-. 4160000000E-01 

0. 

-2.330000000 
1.000000300 
27T7000  080(5 
-.35700e0O00E-01 
49.10000000 

0. _ 

■^755026728566-02 
• 3626869048E-03 
.7912792594 
•8917931248E-01 

“1978434172 - 

-. 133427171 IE-01 
4.419801283 
.2287260839 


. 38800000006-01 


. 1926755760E-02 
. 3756696670E-02 
-. 87 1 631534 IE-02 
♦9997011628 


-4- 

CASE— T-.i 


YF-16  CCV  Canard  FCS  Loops  Closed 


AHAT(ROU  D- 
ArtAT(ROU  2» 
ArtAT ( ROU  3)« 
ArtAT (ROU  4). 

Tflwr(ROu  nv- 

BrtAT(ROU  2)« 
BrtAT ( ROU  3  >• 

JrtATlFOU _ 4»i_ 

FflAT (ROU  !)• 
FflAT  (ROU  2)* 
FflAT  (ROU  3)* 
FflAT  (ROU  4  )• 

"WAT (ROU  I  )•“ 
GflAT  (ROU  2  !• 
CflATfROU  3  )• 
GflAT (ROU  4  )• 


-3.105000000 

-.9796300000 

6.564000000 

0.  _ 

“-37925000000 
. 2988OO0OO0E-O1 
7.178000000 

.0i,.  - 

.6913530975 
-.8143829<7eE-0l 
.5050937166 
.2782O3O729C-01 
“3369331513  “  * 
.20231 12757E-01 
.3508569708 
•2316227340C-91 


8.917000000 

-.2965003000 

-47.71000000 

_0. 

'4.766O0O0e0~ - 

.  18O2O?OOCeE-01 
5.655000000 

-0._ _  _ 

. 7347 1 1 3229 
.9336333501 
-2.79314750? 

- . I  69602 1813 

.3982455362' - 

-.1874449879E-01 

.4955675280 

.2S8S590160E-01 


-.2720000000 

0. 

-8.359000000 

__1 . 00O0O0P00 _ 

1 . 897000000 
- .  3360OO0000E-01 
49.60000000 

JB., _ 

-.1445263004E-01 
.  1026091597E-02 
.356189123? 

_ !_6479790E90E-01 

.1136179162 
- .962S423529E-02 
3.272450037 
.1892645608 


. 3795000000E-0I 
-25.00000000 


.243801 4138E-01 
.2738524359E-02 
-1.626384125 
.9054019678 


NAVION 


CASE  I  T-1/1S 


m 


11- 

-.7770000000 

4.680000000 

21* 

-1.000000000 

-.3556000000 

31- 

1.270000000 

-18.80000000 

4  )• 

6. 

0. 

BflAT(ROU  8 )* 
BdAT (ROU  3)* 

.  BdAT(R0U  4). 
FdAT (ROU  iTT" 
FflAT(ROU  81- 
FdAT (ROU  31- 
FdAT { R0U4) ■ 

“WSTTROU — n~ 
GdAT(ROU  81- 
CHAT (ROU  3). 
CHAT ( ROU  41- 


.78SOO0O0O0E-01 

.7700000000 

e. 


.9383 143219 - 

-.6393779430E-01 
. 9034754 103E-01 
S958S93533E-08 
- .  39465T5033" 

.  1797872883E-01 
.8 17070394 IE-01 
. 1 047484 156E -08 


3.410000000 

.8370000000 

0. 

0. 

“T309S458331 - 

.9661950935 
-.6675638125 
-.841 84467 1 2E-01 
.  1584265886 
. 1034727325E-01 
.1398738707E-08 
.8422088064E-04 


-.4380000000 

0. 

-6.600000000 
1 .000000000 
- .  3T40000W0 

0. 

81.00000000 

0. 

~=T5g53^gg?39E-01 
.  1 13880 3461 E-02 
.6488389070 

_ i.539070699.19^01_ 

- . 37466S8317E-01 
. 12329079S6E-02 
1.131119438 
. 40483S6976E-0  1 


0. 

. 1780000000 

0. 

0. _ _ 


“Tl  7884 1 453T4TT-08 
. 1 129809760E-01 
- .  4159728345E-02 
L99S_9032£A0__ 


CASE  III  --  T-l/15 


YF-16  CCV  Canard  Bare  A/F 


B 


AdAT  C  ROU  11*  -.4310000000 

AdAT (ROU  81*  -1.000000000 

AdAT (ROU  31-  1.670000000 

AdAT  (  ROU  4 )»  0 . 

BdATTROU  IT-  =47170000000 - 

BdATtflOU  81-  .318O000000E-01 

BdAT ( ROU  3  >•  7.630000000 

_  BdAT (ROU  41-  0. 

FMATfRTOU  H-  "70404344356 - 

FdAT(ROU  2  )•  -.6454986763E-01 

FdAT (ROU  31-  .8048464480 

-^gTCROU  4).  .588O856384E-08 

7JMATTR0U IT* =72719177536 
Gd  AT ( ROU  21*  .111 9295553E-01 

GdAT(ROU  31-  .4435147884 

GdATCROW  4)«  . 1554396657E-01 


10.20000000 

-.3060000000 

-50.00000000 

0. 

530000000 - 

. 191O000000E-01 
5.930000000 

0. 

■■76637364488 - 

.9575563303 
-3.995950875 
-.1036480637 
.3012317087 
- . 8788073625E-02 
.3910320437 
. 1302126060E-91 


-.4160000000E-01 

0. 

-2.33000O000 

1.000000000 

8T17 000^000 - 

- . 3570000000E-O1 
49.10000000 
0; _ 

73493658T0 0E=03“ 
.  1671091565E-03 
.8558111402 
_JS174_066395E-01 


.1364919874 
- .6878441879E-02 
3.047796792 
.1040095888 


0. 


. 3880000000E-0 1 


0. 

0. 


.2541 175261E-02 
-.4081544870E-02 
.9999099099 


CASE  II  - 

-  T- 

1/15  YF’16 

CCV  Canard  FCS  Loops 

Closed 

AdAT (ROU 

11* 

-3.105000000 

8.91700000O 

- .8720000000 

9. 

AdAT (ROU 

21* 

-.9796000000 

-.2965000000 

0. 

. 37950O0000E-01 

AdAT (ROU 

31* 

6.564000000 

-47.71000000 

-8.35900C000 

-25.00000000 

_AflATCROU_ 

_4_1*. 

0*.  .  . 

0. 

1 .0O000C000 

0. 

BdAT(ROU 

11* 

-3.985000000 

4.766000000 

1 . 897000C00 

BdAT (ROU 

81* 

.89880000006-01 

.18020000006-01 

- . 3360000A00E-O1 

BdAT (ROU 

3  )• 

7.178000000 

5.655000000 

49.60000000 

^BriAT(ROU 

41* 

-0-. _  _ 

0. 

e. 

FdAT (ROU 

ri¬ 

775354406^5 

.  549"4393'477 

-. I20728T251E=01 

TT3343r98fi6e=?r 

FdAT( ROU 

al* 

-.5792803295E-01 

.9620117449 

. 5219054695E-03 

.21 79316420E-0E 

FdAT (ROU 

31* 

.3710752075 

-2.248227289 

.5314737571 

-1.255750750 

^FdAT(F0U_ 

_4)i 

_ 1 131217993SE-01 

-.846797051  IE-01 

.5010148620E-01 

.9538412096 

CdAT(POU 

1  >• 

-.2373145125 

'  .2825025632  "  ” 

~78965247020r=01 — 

CHAT (ROU 

21- 

.9948800860E-08 

-.8340187072E-08 

-. 5406877 956E- 08 

CHAT (ROU 

31- 

.3055951759 

.3443364719 

8.512771007 

GHAT (ROU 

41- 

. 1199475075E-01 

.  1 181732408E-01 

.9286478638E-01 

75 


NAyiON 


CASE  I  — 

T-. 

94 

flAT ( ROU 

1)- 

-.7770000000 

AHAT(ROU 

2)» 

-1.000000000 

AP1AT(R0U 

3  )• 

1.270000000 

AMAT(ROW 

4  )• 

0. 

“TTOTTffOU  IT:  ^67T5M?M90 
BMAT(ROU  2>»  .72SO0OO00OE-01 

BflAT(ROU  3)»  .7700000000 

BflAT(ROU  4>»  0. 

""FtfATTROU — m - ^9553025756 

FMATtROU  £)•  -.3904748O90E-01 

FriAT(ROU  3)-  .531 1019486E-01 

FMATfROU  4)»  .1048205181 


Gf1AT(ROU 


7683918776E-08 

1999808306E-01 

4712947284E-03 


4.680000000 

-.3556000000 

-12.80000000 

■*2.410000^0" 

.2370000000 

0. 


.98213260S2 
4417088909 


9S671408S2E-01 
.7503300938E-02 
.32S7385388E-03 
.337661 3923E-05 


-.4320000000 

0. 

-6.600000000 

1.000000000 _ 

-.3140000000 

0. 

21.00000000 

=7l49I3lST2«E=5r 
.43756883S8E-03 
7675373788 

4569.32L4Er.01_ 
1 893215507E-01 
.372861491  IE-03 
.7378144217 
. 1540865728E-01 


0. 

.1720000000 

0. 

0. _ 


.682271 1531E-02 
-. 1596940856E-02 


CASE  III  —  T* .04 


YF-16  CCV  Canard  Bare  A/F 


AWATCROU  1  )• 

A7!AT  (RCU  21* 
AP1AT(R0U  3 )• 
.AflAT(ROU  41- 

‘wramrau- tt=" 

BdATlROU  2)- 
BtlATtROU  3 )• 
BDATtROU  4 )■ 

FHAftWCiU  1  >•" 
Ff1AT(R0U  2>« 
FHAKROU  3 )• 
FUAT ( ROU  4)- 
"CTIATTROU — TTT 
GHAT (ROU  2  >• 
GI1AT  (ROU  3)« 
GrtATCRO U  4 1* 


-.4310000000 
-1.0CCC 00009 
1.670000000 

-4. 170W0M0 
•3180000000E-01 
7.630000000 


- • 3930584452E-0 1 
.1013993812 
. 1 80305G655E-02 
"=71649053236 
. 4566014355E-02 
.2826611323 
.5808274716E-02 


10.20000000 

-.3O€0.''0Ot'00 

-50.00000000 

0. _ _ __ _ 

“4.62000SW0 
. 1910300008E-01 
5.930000000 


.9797568870 
-1.879432731 
-.3839724481E-01 
71826802021 
- . 289261 8326E-02 
.2340721474 
. 4694081 842E-02 


-.41 60000000E-01 

0. 

-2.330000000 

1.000000000 

2.170000000 

-.3570000000E-O1 

49.10000000 

0.  _ 

"-TT566478786E-02 
. 61840571 G6E-04 
.9109187967 
.38191 46932E-01_ 
T8'396093475E-01 
-. 30912 13386E-02 
1.880484558 
. 38 1 54434 11E-01 


.  3880000000E-01 


73X591887216-03 
. 1538349980E-02 
.  1 4898 1 3099E-02 
. 9999799274 


YF-16  CCV  Canard  FCS  Loops  Closed 


APIATCROU  1  )•  -3.105000000 

Af1AT(R0U  21*  -.9796^00000 

AflAT<ROU  3  )•  6.564000000 

.  Afl^TtRCU  4  )■  0. _ 


Dim  I  iriWV  •  »  -  --  -- 

8flAT<  ROU  2 )•  .29880000006-01 

BDATIROW  3)»  7.178000000 

_ BrtATCROU  4  )• _ 0;. _ 

^irarrRou  rr=  7B75sr42i84 

FMATCROU  2)*  -.36521230636-01 

FUAT (ROU  3)»  .2389593799 

FMAT(ROU  4 )»  .49386207136-02 

"WTTROU  1  >•  ~  .7*83598848 

CHAT (ROU  2)»  .41427195726-02 

CHAT (ROW  3)«  .2219255433 

CHAT (ROW  4)*  .48581243066-02 


8.917000009 

-.2965000009 

-47.71000009 

~477S600000'0 
.  180200000OE-91 
5.65S000000 

0.  __ _ 

— 73473154457 
.9814276249 
-1.556608858 
-.3339426577E-01 
— H  778481H8 
-.28344508286-92 
.2138328753 
.43556575 l 7E-02 


-.2720000009 

0. 

-8.359000000 

1.000000000 _ 

— 17897000000 
-  .33600000C0E-O1 
49.60000000 
0. 

""=7857771 98706=02" 
.2087635817E-03 
.6985810643 
.3375706817Er01_ 
— T61685773O46-01 
-.2608657566E-02 
1.684841571 
.35621799266-01 


. 379500O000E-01 
-25.00000000 


749S0022584tr=0'2 
.14331967686-02 
-.8451940165 
.  9820997563 _ _ 


76 


CASE  I  —  T-.05 


NAyiON 


AMATCROU 
AMATCROU 
AMATCROU 
AMATCROU 
BTWTTROV 
Bf1AT(ROU 
BMATCROU 
.  BMATCROU 


FMATCROU 
FMATCROU 
FMATCROU 
TKaTCROU 
CMATCROU 
CMATCROU 
CHAT (ROW 


11*  -.7770000000 

2  )•  -1.000000000 

3)*  ». 270O00000 

4>»  e. 

21*  .7250O0O000E-01 

3)-  .7700000000 


4849401 955E-01 
.  669509364 8E-01 
.164B367882E-02 

“2984'446796 - 

. 110G887586E-01 
.2170181841E-01 
. G808G87877E-03 


4.680000000 
-.  3556P00O00 
-12.80000000 
0 » 

274T 0000000” 
.2370000000 
0. 


.9765307325 

-.5322224088 

-.1416026187E-01 

“II931S8269 - 

•8769265325E-02 
•G165845330E-03 
. 798899997GE-05 


-.4320000000 

0. 

-6.G00000000 

1.000000000 

■^73^0000000“ 

0. 

21.00000000 


I  1 a 


.6G71421219E-03 
.7182447813 
.  425759 4997E-01 
""-72540478969E-0r~ 
. 625587 1934E-03 
.8935773G 1 8 
.2357427800E-01 


0. 

.1720000000 

0. 

0. 


.850722S727E-02 
- .243556S041E-02 
.9999581460 


CASE  III  —  T* .05 


YF-16  CCV  Canard  Bare  A/F 


AMATCROU 

A(1AT(R0U 

AMATCROU 

AMATCROU 

BMATCR'OU 

BMATCROU 

BMATCROU 

— BMATXROU_ 
FMATCROU 
FMATCROU 
FMATCROU 
FMATCROU 

“Cmatcrou" 

GMATCR0U 

CMATCROU 

GMATCROU 


1) -  -.4310000000 

2  )•  -1.000000000 
3)-  1 .670000000 

4>«  0. 

2  )■  . 3180000000E-01 

3)*  7.630000000 

ax- _ e, _ 

1  )•  .9660797195 

2) -  -.4887522366E-01 

33*  .1368891187 

43*  .2991324989E-02 

13-  =72053’658134  - - 

23-  . 6721 352 123E-02 

33*  .3457365581 

43*  . 8952949095E-02 


10.20000000 

-.3060000000 

-50.00000000 

0. 

“4762000  0000 
.  191000O000E-01 
5.930000000 


.9722327856 
-2.311075252 
-.59362671 5 4E-01 
“722750 11297 
-.4738986258E-02 
.2925526513 
.7326948138E-02 


-.4160000000E-01 

0. 

-2.330000000 
1 .000000000 
27170000000 
- . 357O000O00E-01 
49.10000000 
0. 

~-T:  r92 5 524494 5-02“ 
.9564 466 181E-04 
.8898722311 
. 4719503489E-01 

7104002 1563 - 

- .  435827 1 650E-02 
2.325886636 
•5919413131E-01 


.  3880000000E-01 


. 1917067339E-02 
-. 230327 1656E-02 
.9999611015 


YF-16  CCy  Canard  FCS  Loops  Closed 


CASE  II  —  T-.es 


AMATCROU 
AMATCROU 
AMATCROU 
„  AMATCROU 
"BMATCROU" 
BMATCROU 
BMATCROU 
_BMATCROU 
FMATCROU 
FMATCROU 
FMATCROU 
FMATCROU 
“CMATCROU" 
GnATCROU 
CMATCROU 
CMATCROU 


-3.105000000 

-.9796000000 

6.564000000 

0. 

“37S2CT0W00 - 

.29880O0000E-01 
7.178000000 
J0_i.  _ 

.84T4S10734 
-.4482256768E-01 
.2910955121 
.7591677878E-02 
“1826709617 
•6051094183E-02 
.2587699R85 
•7267305823E-02 


8.917000000 

-.2965000000 

-47.71000000 

_0.  _ 

4T76600W0O 
.  1802P00000E-01 
5.655000000 

0_.  _ _ 

.4213275315 
.9747742525 
-1.845387524 
-.5043493463E-C1 
72183750549 
- .  4S83226766E-02 
.2638005832 
.6744890998E-02 


-.2720000000 

0. 

-8.359000000 

_ 1 .1.000000000 _ 

1 .89700C000 
- • 33600000006-01 
49.60000000 

_0. _ _ 

-.  1007440058E-01 
.3 1 3636966OE-03 
.6329338519 
.40411 57983E-01 
773236O0726E-01 
- . 35805 30460E -02 
2.020510386 
. 5417478726E-0  1 


. 37950000O0E-01 
-25.00000000 


“7740 12044 2TE-05 
. 1739461895E-02 
-1.012203501 
.9727991840 


NAVI ON 


CASE  I  —  T-.02 


AHAT(ROU 

1)- 

-.7770000000 

(.680000000 

-.4320000000 

0. 

AHAT<ROU 

2  )• 

-1.000000000 

-.3556000000 

e. 

.1720000000 

<  AflAKROU 

3)» 

1.270000000 

-12.80000000 

-6.600000000 

0. 

£5$J(R0U 

4  )• 

0. 

0. 

1.000000000 

0. 

warmou 

1  )• 

- 6. 10000000 0 

2.(10000000 

-.31(0000000 

BHAKROU 

2)« 

. 725000O000E-01 

.2370000000 

e. 

BHATtROU 

3). 

.7700000000 

0. 

21.00000000 

_  BHAT(ROU 

4)- 

e. 

0. 

e. 

fhakRou 

1  )• 

•9835(51199 

• 9356^3 4665E-0X 

=7S@2BX7(X(5E=52 

7161597  4  1  (72E-03 

FHAKROU 

2  )• 

-. 1976761 620E-01 

.9919767251 

.11(92788506-03 

.342S722506E-02 

FHAT(ROU 

3)- 

.2603042592E-01 

-.2377387318 

•676231(076 

-.(191740962E-03 

•FHAKROU 

4). 

.2582943506E-03 

- • 2(37058699 E-02 

. 1873548183E-01 

"CTfSTTROff- 

IK 

“=7X2X0102(58 - 

T(803255323E-0i 

-7795693981  IE-02 

CHAT(ROU 

2)« 

.2655609736E-02 

.42(3(150676-02 

•7860836784E-04 

CHAT (ROU 

3)« 

. 1267(038726-01 

.4(9064 73 3(E -04 

.39336(01(1 

GHAT (ROU 

4  )• 

• 13S7845151E-03 

. 2350923 155E-06 

. (020525562E-02 

YF-16  CCV  Canard 

Bare  A/F 

CASE  III 

— ■  T« 

.02 

AHAT(ROU 

U* 

-.4310000000 

10.20000000 

-.41600000006-01 

0. 

AHAKROU 

21* 

-1.000000000 

-.3060003000 

0. 

. 3880000000E-01 

AHAKROU 

3)- 

1.670000000 

-50.00000000 

-2.330000000 

0* 

AflAT  ( ROU 

4  )• 

0. 

0. 

1.000000000 

0. 

BHATIROU 

1  !• 

-4.170003000 

~ — 47620W0000 - 

2. 170000000 

BHAKROU 

21* 

.31808000006-01 

.  1910000000E-01 

- . 3570000000E-OI 

BHATtROU 

31« 

7.630000000 

5.930000000 

49.10000000 

—tMATlROU. 

4). 

0. 

0. 

0. 

fWaTTrou 

1  >• 

.98937724(3 

.  2027736364 

80831 10994E-03 

.788(31598(6-04 

FHAT(ROU 

2  >• 

-.19839(47186-01 

.9918682(80 

. 1577175559E-04 

.77310(42286-03 

FflAKROU 

31- 

•4226(537(CE-01 

-.97003972(1 

.95(450(6(7 

-.38026255966-03 

_ FflAT  (RCU_ 

4  )• 

.  ,393(5307256-03 

-.9800581(326-02 

.195(1038(?E-@1 _ 

.999997(620 

CHAT (ROU 

1  )• 

-.829824361 IE-01 

.9193036(136-01 

.(271056099E-01 

GHAT (ROU 

2)- 

.1(6(0607226-02 

-. 53795(71 (gc-03 

-.1137872572E-02 

GHAT (ROU 

3  >• 

.1(71(57657 

.1175089202 

.9606686627 

GHAT'kOU 

4>» 

.1489947S93E-02 

. 1178197980E-02 

.9677019529E-02 

YF-16  CCV  Canard  FCS  Loops  Closed 


CASE  II  —  T..02 


AWAT  (  ROU 
APIAT  <  ROU 
ArtATCROU 
ADAIiROU 
BriAT  ( ROtr 
BHAKROU 
BOAT (ROU 
«HmT(ROU 
FflAT  ( ROU~ 
FHAKROU 
FHAKROU 
— E.HAT  (  ROU 
’  CflAKROlT 
CHAT ( ROU 
CHAKROU 
GHAT (ROU 


1  )■  -3. 105000000 

2) «  -.9796000000 

3) *  6.56(000000 

4  >•  0. 

Tl-  =37925000000 - 

2>«  .29880000006-01 

3) -  7.178000000 

4) *  0, 


H-  :937777028s - 

2  )•  -. 189256 81 68E-01 

3) -  .125(679558 

4) «  . 127(224(136-02 

i  )•  -77636796768E-01- 
2>»  .13(978784  IE-02 

3)*  .1266971358 

4  )•  .  1322227582E-02 


8.917000000 

-.2965000000 

-47.71000000 

0. 

47T6E00?g00 
. 1802000000E-01 
5.655000000 
0. 

7X7(6757(79 - 

.9923652712 
-.8631235396 
-. 89295 61842E-05 
.921O287677E-01- 
-.5512878322E-03 
.109877(464 
. 1109280185E-02 


-.2720000000 

0. 

-8.359000000 
_ L.  0000000.00_ 


1.897000000 

-.3360000000E-O1 

49.60000000 


0. 


rX8‘40'832552T-'0e 
.5640188019E-04 
.8(12450672 
.  18384G8656E-01 
73(203015816-01“ 
-.  1013843279E-02 
.91(5976885 
.93995 19380E-02 


0. 

.3795000000E-01 
-25.00000000 
JU _ 


TT325ff9(92rr=B2_ 
.7(67331 653E-03 
.4599560387 
.9952692512 


78 


NAVION 


CASE  I  --  T».®125 


By 


h 


AHAT(ROU  i)« 
AHAT(ROU  2)- 
AHAT(ROU  3)« 

_ _ flftAT(ROU  4). 

BHATTRCU — mr 
BRAT ( ROU  2 ) ■ 
BHAT(ROU  3)- 
_  8HAT(ROU  4  )• 

rrrsTTRog — fv 

FHAT(ROU  2)- 
FHAT(ROU  3)» 
FUAT (ROU  4>- 
~criAT(Rog — ttf- 
GHAT(ROU  2)- 
CHAT (ROU  3  )• 
CHAT (ROU  4>- 


-.7770000000 
-1.000000000 
1.270000000 
0. 


-. 124 1007 876E-01 
.  1G1 2929741E-01 
.1O0292S153E-03 
"^77586926797EH0r 
. 1378562340E-02 
•855S443942E-02 
• 55697464 12E-04 


4.680000000 
-.3556000000 
-12.80000000 
0. 


•7250000000E-01  .2370000000 

.7700000000  0. 

0.  0. 


r5^rg7^ME^0r 

.9951990693 
-.1527613544 
9696524 003E -03 
73006 f75T47EH0'l 
•  2768182S1  IE-02 
.U89734304E-04 
. 3953238250E-07 


-.4320000000 

0. 

-6.600000000 
1.000000000 
"=731T0’000-000- 

0. 

21.00000000 

0. 

I4573793149E-04 
.9207685898 

_ sJ  1 99808842E-01 

- .  459281 69552^02“ 
.2844749554E-04 
.2519283650 
.1 5962806 16E-02 


0. 

.1720000000 

0. 

0. 


•2144966566E-02 
-. 1667802128E-03 
.9999992997 


CASE  III  —  T-.012S 


YF-16  CCV  Canard  Bare  A/F 


IN 

,'4 


AH  AT  ( ROU  1  )- 
ARAT ( ROU  2)> 
AHAT(R0U  3 >• 
%iAnAT  ( ROU  4  1« 


bhat ( Rou  rr 

BRAT  ( ROU  2>- 
BHAT ( ROU  3>- 
MATJ.ROU  4 )- 

tMAT  1  DAI  I  rr~ 


-.4310000000 

-1.000000000 

1.670000000 

0. 


’THAT (ROU  n* 
FUAT <  ROU  2)* 
FHAT(R0U  3)- 
_£TIAT(R0U  4). 
CRATTROU  1  )•" 
CHAT (ROU  2)0 
CHAT (ROU  3 )■ 
CHAT (ROU  4 )• 


"--C17000W00 - 

.3180000000E-O1 

7.630000000 

0. 


10.20000000 
-.3060000000 
-SO. 00000000 

0. 


.9938280622 
-. 1243920345E-01 
.2436901531E-01 
.  1 4508 1 3907E-03 
- . 51 9701 7096E-0 1 
•7215721478E-03 
. 93271 62452E-01 
. 5873774656E-03 


4.620000000 
. 1910000000E-01 
5.930000000 
JU. 


.1270405060 
.9953876977 
-.6133245727 
VJ.3857651595E-02 
.  57575  45600E-0r 
•.  121398 1239E-03 
. 736523974 IE-01 
.46 1230636 IE-03 


- . 4160000000E-01 

0. 

-2.330000000 

__1 . 000000000 _ 

2TT70000000 
- . 3570000000E-01 
49.10000000 
-0. 


0. 

. 3880000000E-01 

0. 

0. _ 


- .  51 0850431 8E-03 
.6206153958E-05 
.9712884452 
.  1231969689E-01 
.2685873820E-01- 
-.6129946312E-03 
.6053496619 
. 3800817208E-02 


CASE  II  —  T-.0125 


YF-16  CCV  Canard  FCS  Loops  Closed 


AHAT(ROU  1>« 
ARAT (ROU  2)« 
Af1AT(R0U  31* 
-AHAT.IROU  4). 
' BRAT (ROU O' 
BRAT (ROU  2)« 
BRAT (ROU  3)« 

_ BRAT(ROU  4  )• 

FRAT(ROU 
FUAT (ROU  2 )« 
Ff1AT(R0U  3>- 
FRAT(ROU  4)4 
— CHAT  ( ROU — n< 
CHAT (ROU  2 )« 
CHAT (ROU  3>* 
CHAT (ROU  4 )• 


-3.105000000 
-.9796000000 
6.S64000000 
0. 


“37925000000 
.2988000000E-01 
7.178000000 
0. 

-!l 198467995E-01 
.7979469095E-01 
•S034431973E-03 
“48234156135-01 
. 669538 3975E -03 
.8305417439E-01 
.  5328398624E-03 


8.917000000 

-.2965000000 

-47.71000000 

0. 

"”47766000000" 

.  1802000000E-1 
S.  655000000 


.9956241283 
-.5603306928 
- . 357634010  IE-02 
158315615605-01 " 
-. 1 342S53646E-03 
.6940796792E-01 
.4364282181E-03 


-.2720000000 

e. 

-8.359000000 
1 .000000O00_ 


57897000000 
01  -.3360000000E-01 

49.60000000 

0. 


■^T5l62405534E-02 
. 2267550034E-04 
.8988316140 
.1186083171E-01_ 
72222406262E-01~ 
-.5575070930E-03 
.5893724494 
.  3746696550E-02 


.3084586472E-04 
.4839453630E—03 
- . 14967688 19E-03 
.9999993744 


0. 

. 37950OOO00E-O1 
-25.00000000 
As. - 


. 532562381 6E-03 
. 4709993003E-03 
-.2966565147  . 

_ « 998.LL32422 - 


79 


NAVION 


CASE  I  —  T-.Ol 


ARATCROU  1  )• 
ARATCROU  2  )• 
AflAT < ROW  3 )• 
LOW _ 4  I  • 


-mm 


. . row  ri=" 

SHAT (ROW  2)* 
BRAT (ROW  3)* 

-ffifBS-fR- 

FUAT (ROW  3)« 
FUAT (ROW  3>« 

. FflATtROU  4  >■ 
CHAT  ( ROd — TTS" 
GHAT (ROW  2)- 
CHAT (ROW  3)* 
CHAT (ROW  4)> 


-.7770000000 

-i.eoeopeooo 

1.870000009 

-rrmewaoo - 

.725OOOOOOOE-01 

.7700000000 

JL, 


4.680900000 

-.3556000000 

-12.80000000 


TS9S0MM36 - 

-.9942633432E-02 
. 1286448839E-01 
. 64054 13975E-04 
r76075?59'49!rrO'r 
. 1027569263E-03 
. 70153291 78E-02 
. 362151144SE-04 


.2370000000 

0. 

0. 


T«FOTi?T7re=0r 

.9962162200 
-.1233444327 
-.6244045283E-03 
24O5998227E-0T 
. 22455655 80E-02 
.6386603604E-O5 
.171 4434361 E-07 


-.4320000000 

0. 

-6.600000000 

_ L..O.00OO00O0_ 

-.3140000000 

0. 

21.00000000 

0. 


0. 

.1 

0. 

.A. 


.1720000000 


.2945506434E-04 

.9361034362 

_ j967705033?E-O2 

~-73570rS330’6rr02' 
•  1771492273E-04 
.2031979441 
. 1027204472E-02 


. 4025012640E-04 
. 171681 1258E-02 
-. 1073975789E-03 
.9999996398 


YF-16  CCV  Canard  Bare  A/F 


CASE  III 

• 

►- 

I 

« 

.01 

ARATCROU 

1>- 

-.4310009000 

10.20000000 

- . 4160O00O00E-01 

0. 

ARAT ( ROW 

a>- 

-1.000000000 

-.3060000000 

0. 

. 388000OOO0E-O1 

AflAT  ( ROW 

31- 

1.670000000 

-50.00000000 

-2.33000eO0O 

0. 

4  )• 

0. 

0. 

1 .000000000 

0. 

r~r  ii  i  r  r  n  — — 

1 1 1 1  i  1 1  — 

bhatcrou 

2)« 

. 3180000000E-01 

. 191O0000OOE-O1 

-. 3570000O00E-01 

BRAT (ROW 

31“ 

7.630000000 

5.930000000 

49.ieoeeeoe 

BOAT (ROW 

4)- 

0. 

0. 

0. 

TV*  ‘ 

■  1 1 1 II II  — 

mm  i  mi  1 1  — 

KKiUltimiCKH 

FRAT  <  ROW 

2)« 

-.99615009456-02 

.9964358261 

.3981578242E-05 

.387341 U56E-03 

FRATCROU 

31- 

. 1894283792E-01 

-.4925337505 

.9769652730 

-.960351 51 82E-04 

FHATCROU 

4)- 

. 9099833036E-04 

-.247S132779E-02 

. 9884386482E-02 

.9999996791 

CHAT (ROW 

!)• 

-.4160268085E-01 

•460901974SE-01 

.2153024875E-01 

CHAT (ROW 

21* 

.5255320796E-03 

- . 396S890725E-04 

- . 4639650323E-03 

CHAT (ROW 

3). 

. 749596 9660E-O1 

•5898754909E-01 

.4856092049 

GHATtROU 

4)- 

. 37705 1 9690E-03 

.295425025 IE-03 

. 243698 1830E-02 

YF— 16  CCV  Canard  PCS  Loops  Closed 


CASE  II  —  T-.Ol 


ARATCROU 
ARAT ( ROW 
AflATfROU 
AHAHROW 
BTI  AT  (ROW 
BRAT (ROW 
BRATCROU 
BOAT (ROW 
FWaTCROW 
FRATCROU 
FRATCROU 
FRATCROU 
CHAT(ROU~ 
GHAT(ROU 
CHAT (ROW 
CHAT (ROW 


D-  -3.105000000 

2) «  -.9796000000 

3) «  6.564000000 

41-  0. 

rr-  -3 7925000000 - 

2) *  . 2988O0O0O0E-O 1 

3) -  7.178000000 

4) «  0. 


8.917000000 
-.2965000000 
-47.71000000 
0. 


T7S 79689723375 

SJ*  - . 9629390870E-02 
3)»  .6419997708E-01 

J >i___i32341 1S495E-03 
ITS -7387208460OE-01~ 

2) «  .4887348029E-03 

3) »  .674927 171 5E-01 
4>-  . 3445490309E- 03 


47766009000 
.  1802009009E-01 
5.655000000 
_0U_ 


.  8827758866E-01 
.9966056658 
-.4539311793 
_r_.23O7934l66E-02_ 
.  4685342750E-01 
-.5060593333E-04 
.55726551 78E-01 
•2799897326E-03 


-.2720000000 

0. 

-8.359O00OOO 

_JL.00090O000_ 


1.897000000 
-.  3360O00000E-O1 
49.60000000 
0^ 


-.25671 76464E-02 
. 1465 1 71999E-04 
.9185402417 

_ ,3589 1 70877ErO2_ 

. 1801 119674E-01 
- .  4248695386E-03 
.4763139338 
•2414O96702E-O2 


0. 

. 3795OO090OE-O 1 
-25.00000000 


.344O086304E-O3 
.3776S02744E-03 
-.2398168580 
_ .9987840890 _ 


80 


APPENDIX  B 


VBA  HAnnWABK  AND  SOFTWARE  SYSTEMS 

VBA  SYSTEMS  OVERVIEW 

The  Princeton  Flight  Research  Laboratory's  Variable  Response 
Research  Aircraft  (VRA)  (Figure  24)  consists  of  a  modified  NAVION  air¬ 
craft  with  independent  control  of  lift,  sideforce,  thrust,  pitch,  roll 
and  yaw.  Mediating  between  the  evaluation  pilot  and  actual  aircraft 
control  are  a  specially  constructed  analog  computer  system  and  a  digital 
microprocessor  (Micro-DFCS) ,  which,  in  the  present  research  application, 
accepts  8  channels  of  analog  sensor  and  command  data  and  supplies  three 
channels  of  output  control  which,  in  turn,  drive  control  surfaces  actua¬ 
tors. 

Independent  control  of  three  forces  and  three  moments  is  provided  by 
commands  to  the  elevator,  ailerons,  rudder,  throttle,  direct-lift  flaps 
and  side-force  panels.  The  control  surfaces  are  driven  by  hydraulic 
servos  capable  of  high  surface  accelerations  (some  were  originally 
fitted  to  the  B-58  aircraft).  The  modified  VRA  units  incorporate  sole¬ 
noid  actuated  valves  with  force-override  features  for  quick  disen¬ 
gagement.  Characteristics  of  the  control  surface  effectors  are  sum¬ 
marized  in  Table  6.  Surface  rate  limits  are  seen  to  range  from  60  to 

110  deg/sec.  Bandwldths  are  given  for  flat  response  and  6  db  attenua¬ 
tion  (in  parentheses),  except  that  the  thrust  bandwidth  is  specified  by 
the  frequency  for  3  dB  attenuation.  The  aircraft's  normal  operating 
speed  range  is  65  to  120  kt;  maximum  specific  forces  and  moments  ("con¬ 
trol  power")  are  given  for  70  kt  airspeed.  At  IAS  «  105  kt  maximum 
direct  lift  and  side-force  accelerations  are  1  g  and  0.5  g,  respec¬ 
tively. 

SENSORS  AND  COMMAND  SIGNALS 

The  sensors  used  for  most  flight  testing  include  angular  rate  gyros 
and  linear  accelerometers  for  all  three  axes,  vertical  and  heading 
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gyros,  dual  angle-of-attack  and  sideslip  angle  vanes,  radar  altimeter, 
indicated  airspeed,  control  surface  positions  and  cockpit  control  posi¬ 
tions. 

In  conventional  application  of  the  VRA  (Figure  25a),  aircraft  sensor 
signals  and  command  signals  from  the  evaluation  pilot's  controls  are  fed 
into  the  VRA  analog  system  which  can  be  set  with  the  necessary  gains  and 
crossfeeds  to  emulate  performance  characteristics  of  a  variety  of  other 
aircraft.  The  .licro-DFCS  was  originally  conceived  as  a  research  tool 
for  the  evaluation  of  digital  systems  in  flight  control.  As  presently 
implemented  in  the  VRA,  the  Micro-DFCS  operates  in  parallel  with  the 
analog  system.  Although  capable  of  exerting  full  6-axis  control,  the 
digital  system  is  usually  enabled  in  only  3  axes,  with  the  analog  system 
supplying  the  remaining  vectors. 

For  the  present  research  application,  the  VRA's  analog  simulation 
and  the  Micro-DFSC  serve  complementary  roles.  The  analog  system  is  set 
to  emulate  the  lateral  directional  characteristics  of  the  bare-airf rame 
YF-16;  however,  the  gains  mediating  pilot  commands  through  it  are 
reduced  to  zero.  Thus,  all  of  the  command  augmentation  for  the  present 
problem  is  supplied  by  the  digital  system,  while  the  analog  system, 
responding  only  to  feedback,  maintains  the  bare-airf rame  YF-16  configu¬ 
ration.  This  relationship  is  summarized  graphically  in  Figure  25b. 

MICROPROCESSOR  DIGITAL  FLIGHT  CONTROL  SYSTEM 

The  Micro-DFCS  has  undergone  substantial  evolution  since  the  begin¬ 
ning  of  the  contract  period,  due  largely  to  FRL's  recognition  of  defi¬ 
ciencies  in  the  original  system  in  terms  of  throughput  delay  and  program 
memory.  They  have  also  expanded  program  development  capability  and 

facilities.  Concomitantly,  STI  has  become  increasingly  sophisticated  in 
the  use  of  the  software  development  resources  available. 

Initially,  the  Micro-DFCS  was  configured  with  Intel  SBC  family 
boards.  The  processor  board  (iSBC  80-05)  contained  an  8085  microproces¬ 
sor  running  from  a  2  mHz  system  clock.  In  order  to  perfoim  floating 
point  operations,  an  Intel  iSBC-310  ,?high  speed  math  unit"  (a  separate 
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a)  VRA  System  Modules 


b)  DFCS  Module  Interconnections 

Figure  25.  VRA  Analog  and  Digital  Systems  —  Data  Path  Diagrams 
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multibus  board)  was  also  Installed.  In  addition,  several  boards  of  RAM 
(memory)  were  also  required.  Control  software,  a  6K  byte  assembly 
language  program  (CAS-1)  was  written  to  operate  in  this  environment. 

In  September  1980,  the  Intel  80-05  board  was  replaced  with  a  Mono¬ 
lithic  Systems  80-04  processor  board  containing  a  Z-80  microprocessor 
with  32K  bytes  of  onboard  RAM,  a  4  mHz  clock  and  a  single-chip  floating 
point  math  unit  (AMD  9511),  with  capabilities  considerably  beyond  the 
iSBC-310  math  unit.  By  this  time,  the  CAS-1  control  program  had  evolved 
to  CAS-4  which,  like  its  predecessors,  still  accessed  the  iSBC-310  math 
board  for  floating  point  calculations.  An  upgraded  version  called  CAS-6 
was  written  for  the  new  Z-80  based  system;  however,  the  enormous  incon¬ 
venience  of  adapting  the  assembly  code  to  accommodate  the  different  cod¬ 
ing  convention  required  by  the  AMD-9511  math  chip  dictated  that  CAS-6 
also  use  the  ponderously  slow  iSBC-310  math  board. 

SOFTWARE  DEVELOPMENT  SYSTEM 

STI's  initial  adaptations  of  FRL's  software  to  the  YF-16  problem 
centered  around  the  CAS-6  program.  However,  when  certain  subtle  dis¬ 
crepancies  in  the  responses  of  the  YF-16  zero-order-hold  problem  were 
observed,  an  extensive  static  check  of  the  longitudinal  modes  of  CAS-6 
was  performed  by  J.  Smith  and  T.  Myers  of  STI.  In  the  course  of  this 
check,  it  was  discovered  that  a  serious  "scramble"  existed  in  the  gain 
selection  algorithm  of  CAS-6  and  that,  apparently,  the  program  had  not 
been  thoroughly  static-checked  by  its  original  author.  This  discovery 
hastened  FRL's  decision  to  abandon  the  assembly  language  code  and  pursue 
the  installation  of  higher  level  language  software  which  would  support 
the  onboard  floating  point  capability  of  the  MSC  8004  hardware. 

Up  to  this  point,  program  development  for  the  CAS  series  of  assembly 
language  programs  had  been  accomplished  using  a  FORTRAN  cross  assembler 
resident  on  Princeton's  IBM  3033  timesharing  system.  Assembly  code 
would  be  written  using  the  3033  system  editor  facilities  and  the  result¬ 
ing  source  code  assigned  as  an  input  file  to  one  of  several  disk  resi¬ 
dent  cross  assemblers.  The  cross  assembler,  in  turn,  created  an  output 
file  in  standard  Intel  hex  format  which  was  then  loaded  on  to  the  80-04 


microprocessor  over  a  standard  voice-grade  phone  link  at  300  baud.  When 
the  relative  Inefficiency  of  this  system  became  obvious,  FRL,  with  some 
counsel  from  ST1,  made  the  decision  to  Implement  a  local  disk-based 
development  system  using  a  processor  compatible  with  the  Z-80  on  the 
80-04  board.  They  acquired  a  MSC-80-09  processor  board  from  Monolithic 
systems,  which  represented  one  of  the  most  advanced  single-board  compu¬ 
ters  available.  The  80-09  contained  a  Z-80  processor,  9511  math  unit, 
floppy  disk  controller,  ROM  with  bootstrap  for  a  standard  CP/M  operating 
system,  two  serial  RS-232  ports  and  32K  of  RAM  on  a  single  multibus 
board.  Additional  elements  included  dual  8  in.  floppies  (adapted  from 
an  earlier  system),  a  LSI  ADM-31  CRT  terminal,  and  an  Anadex  DP-9501 
matrix  printer. 

MWITOR  AND  LANGUAGE  SELECTION 

The  selection  of  CP/M  as  the  operating  monitor  for  the  development 
system  opened  a  choice  among  a  number  of  higher  level  languages  which 
could  be  used  to  c^nerate  Z-80  code  acceptable  to  the  MSC  80-04  board 
which  remained  as  the  heart  of  the  Micro-DFCS.  Two  versions  of  the 
PASCAL  language,  FORTRAN  and  a  compiler  version  of  BASIC,  were  evaluated 
for  suitability  and  benchmarked  for  object  code  efficiency.  The  defini¬ 
tive  winner  was  PASCAL /MT  v.  3.2,  created  by  MT  Microsystems  of  Carls¬ 
bad,  CA.  At  the  time  of  the  benchmark  runs,  PASCAL  /MT  3.2  had  already 
been  obsoleted  by  PASCAL /Mr  v.  5.x.  The  latter,  however,  while  contain¬ 
ing  more  convenience  features  than  its  predecessor,  generated  markedly 
less  compact  object  code.  Even  using  the  less  replete  of  the  two 
PASCAL,  the  execution  module  for  a  PASCAL  equivalent  of  CAS-6  still  con¬ 
tained  2-3  times  the  number  of  instructions  of  the  assembly  code  CAS-6. 
Nevertheless,  the  introduction  of  the  PASCAL-based  system  rendered 
floating  point  operations  and  program  logic  much  more  accessible  than 
was  possible  with  the  assembly  language  code.  This  permitted  much  more 
flexibility  in  terms  of  control  algorithms,  and,  in  general,  saved  con¬ 
siderable  time  and  effort  in  the  development  of  the  revised  ZOH  and 
slower  code. 
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S0FTK4RE  DEVELOPMENT  FOR  SLEWER  CONTROL  LAN 

During  the  course  of  the  contract  period,  two  major  efforts  at 
Slewer  Implementation  were  mounted*  The  first  effort  (February  1981) 
centered  around  a  formulation  which  had  the  capability  of  handling  three 
model  input  vectors.  This  formulation  was  coded  into  an  adaptation  of 
the  FRL  CAS-6  assembly  language  code.  Although  it  was  immune  from  the 
gain  selection  glitch  inherent  in  the  Zero  Order  Hold,  it  nevertheless 
proved  extremely  resistant  to  static  check.  This  served  as  a  beginning 
object  lesson  in  the  validation  of  slewers,  i.e.,  that  careful  formal¬ 
ized  validation  procedures  whose  numerical  performance  can  be  compared 
to  theoretical  check  cases  are  a  prime  necessity. 

The  second  major  onslaught  against  the  Slewer  occurred  after  the 
decision  had  been  reached  to  abandon  the  assembly  language  version  of 
CAS-6  in  favor  of  the  framework  of  R.V.  Walters'  PASCAL  language  code 
(pCAS).  Groundwork  was  laid  by  first  adapting  the  Zero  Order  Hold 
within  the  pCAS  framework,  and  performing  extensive  validation  work  to 
assure  ourselves  that  the  software  was  indeed  performing  to  specifica¬ 
tion.  The  first  PASCAL  slewer  embodied  the  February  1981  Slewer 
algorithm;  however,  that  effort  was  truncated  when  the  excessive  compu¬ 
tational  load  of  that  model  became  apparent.  A  second  (October  1981) 
slewer  formulation  proved  to  be  computationally  more  congenial  for  two 
reasons.  First,  it  utilized  only  pedal  command  input;  second,  it  had  a 
more  balanced  distribution  of  foreground  and  background  computation  than 
its  predecessor.  Rough-out  code  for  this  latter  version  of  the  Slewer 
were  completed  and  a  preliminary  compiled  version  generated. 

Overall  structure  of  the  software  is  shown  in  Figure  26.  The  slewer 
control  program  SLEW2  contains  two  functional  modes  of  operation  which 
can  be  selected  by  the  evaluation  pilot  using  the  cockpit  Termlflex  con¬ 
trol  console.  Mode  "E"  initiates  ZOH  operation  and  mode  "S"  initiates 
slewer  operation.  Because  of  the  RAM  space  required  by  the  additional 
gains  and  constants  associated  with  the  Slewer,  it  was  necessary  to 
strip  out  all  of  the  longitudinal  axis  setup  and  control  routines  as 
well  as  some  of  the  analog  test  utilities  from  the  currently  operational 
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DFCS  Hardware  MSC  80-09  CP/M  System 


Figure  26.  DFCS  Hardware  and  Software 


version  of  pCAS  (ver.  4.3).  In  addition,  several  alternative  versions 
of  the  slewer  control  program  had  to  be  compiled  In  order  to  address 
certain  validation  problems.  These  will  be  discussed  later  on.  It 
should  be  noted  that  the  pCAS  control  program  containing  the  ZOH  con¬ 
troller  also  contained  an  operational  slewer.  However,  the  scaling  and 
graceful  entry  problems  discussed  earlier  In  the  text  made  it  too 
hazardous  to  engage  mode  "S"  during  any  of  our  last  three  flights. 

NUMERICAL  VALIDATION  OF  THE  SLEWER  ALGORITHM 

Due  to  the  dynamic  structure  of  the  Slewer  algorithm,  a  static  end- 
to-end  I/O  check  was  not  possible  using  the  procedures  developed  for  the 
ZOH  case.  It  was  necessary,  therefore,  to  compile  a  special  validation 
version  of  the  Slewer  in  order  to  verify  the  internal  calculations  of 
the  real-time  flight  controller  program  (SLEW2).  The  validation  program 
(SLTST)  differed  from  the  real-time  controller  in  three  regards: 

1.  The  analog  drivers  were  disabled  so  that  the  pro¬ 
gram  could  function  using  the  80-04  processor  in 
a  stand-alone  environment. 

2.  LPRINT  statements  (PASCAL  printer  output  com¬ 
mands)  were  inserted  into  the  control  loop  sec¬ 
tion  of  the  source  code  in  order  to  print  out 
intermediate  and  output  variables  at  every  sample 
frame. 

3.  Unity  values  were  substituted  for  the  front-  and 
back-end  conversion  scale  factors  so  that  values 
in  pure  physical  units  would  be  displayed. 

In  order  to  validate  the  numerical  output  of  the  Micro-DFCS,  a  vali¬ 
dation  case  consisting  of  output  values  for  three  controllers  were 
worked  out  using  the  discrete  model  with  constant  pedal  input  and  fixed 
(constant)  feedback  variables.  Numerical  validation  of  the  October  1981 
algorithm  was  obtained  after  resolving  some  tricky  problems  resulting 
from  uninitialized  variables  in  the  PASCAL  code. 
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ANNOTATED  TOTAL  SDN 


This  run  checks  verisimilitude  of  pseudo  Inverse,  ZOH  matching  of 
continuous  system  by  comparing  s-plane  matrices  (continuous  model) 
against  w'-matrices  (simulator). 
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APPENDIX  D 

VALIDATION  OIF  ZOH  CONTROL  LAN,  FREQUENCY  RESPONSE 

A  comparison  of  the  theoretical  controller  working  against  an  Ideal 
model  of  the  open- loop  aircraft  was  carried  through  using  the  frequency 
response  concepts  of  Reference  2.  A  comprehensive  set,  at  10  Hz  and 
50  Hz,  follows. 


Y-/6  CCV  Lower  Order  Equivalent  System  Model 
Figure  29.  Bode  Plot  10  Hz,  50  Hz 


Figure  30,  Bode  Plot  8/5. 


YF-/6  CCZ  Lower  Order  Equivalent  System  Mode / 


Figure  32.  Bode  Plot  0/ 


YE-/6  CCV  Lower  Order  Equivalent  System  Model 

Figure  33.  Bode  Plot  r/6c 


Figure  34.  Bode  Plot  8/ 6 


YF-16  CCV  Lower  Order  Equivalent  System  Model 


Figure  35.  Bode  Plot  p / Sq 


APPBSDIX  B 


SLEW*  FIRST-OIDBR  FORK 


Given 


x  ■  Ax  +  Bu 


(122) 


the  Laplace  transform  is 

X (s)  -  [la  -  Aj^xCO)  +  [la  -  A]"1  BU  (123) 

If  the  continuous  control  vector  Is  a  reconstructed  signal,  using  slewer 
data  holds,  then 

U(s)  -  M(s)RT(s)  (124) 

where  M  Is  the  a-plane  transform  of  the  data  hold.  For  the  slewer.  If 
T 

R  la  the  sampled  Input, 

-  M1(s)RT  +  M2 (s) (RTe”sT)  (125) 


The  motivation  for  partitioning  Equation  125  in  this  particular 
manner  is  clear  when  the  Impulse  responses  of  and  M2  are  plotted. 


From  Equation  125  write  the  Impulse  responses: 


Mi<t) 


M2(t) 


1  -  (t/T) 

0 


0  <  t  <  T 

Elsewhere 

0  <  t  <  T 

Elsewhere 


(126) 


(127) 


From  Equation  125, 

X(s)  -  [is  -  a]-1x(0)  +  [is  -  A]"1B[M1RT  +  M2RTe“sT] 


Let 


[is  -  A]-1  -  $(t) 


(128) 


(129) 


and  write  the  Inverse  of  Equation  123: 


(/: 


x(t)  -  ♦(t)x(O)  +  I  I  4>(t  -  5)BMi(5)  d5 

*o 


rT(t) 


[/ 


+  11  *(t  -  5)BM2(5)  dC 

o 


rT(t  -  7) 


(130) 


Using  the  definition  of  Equation  126  and  127,  substituted  Into  Equa¬ 
tion  E-9,  yields  the  explicit  solution 


x(t) 


♦(t)x(0)  + 


4>(t  -  5)  $  5  d5 


RT(t) 


+ 


RT(t  -  T) 

(131) 
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Next,  turn  the  Integral  equation  into  an  integro-recurslon  equation  by 
setting 


t-kT  +  x  ,  k  -  0,  1, 


Clearly,  in  the  recursion 


,  0  <  T  <  T 


(132) 


x(2T) 


x(kT)  =£>  x[  (k  +  1)T] 


Substituting  Equation  132  Into  131  gives  the  first-order  form 


/•0  p 

x(kT  +  t)  -  ^(t)x(kT)  +  /  *(t  -  5)  f  C  d£  Rt (kt  +  t) 

L-'O 


*(t  -  5)B[1 


-  U/T)l  dejtB 


RT[kT]  (133) 


Finally,  suppress  the  redundant  notation  with  regard  to  T,  and  thus 
verify  the  slewer  first-order  form 


xk+l  "  $xk  +  rlrk  +  r2rk-l 


(134) 


where 


+  -  [is  -  A]' 


•  i  f 

•'A 


4>(t  -  C)B  d? 


+  -1  !  [Is  -  A] -1b 

T  s2 

s  Jt-T 


(135) 


(136) 


(137 


r2 


1 

T 


t  -  5)B[1  -  (5/T)] 


t-T 


In  going  from  Equation  133  to  Equation  134  we  make  use  of  the  dis 
T 

Crete  nature  of  R  ,  specifically, 


RT[(k  +  1)T]  =  R(kT) 

RT[(k  -  1)T]  =  R(kT) 


(138) 


since  the  sampling  value  is  determined  at  the  start  of  the  frame  time. 


Finally,  observe 


* 


'-«■  r«u«:a. 


If9 


r<? 


ly 

ii 


APPENDIX  F 


SLEWER  VIENSR-HOPF  EQUATIONS 


Given  a  discrete  set  of  equations  of  motion  for  the  model 


5i 

S 


i 


x®k+l  *  +  ^m^mjt 


(HO) 


one  desires  to  modify  the  simulator,  modeled  by  the  slewer  data  hold 
equations 


Xsk+1  "  *sxsk  +  rl6sk  +  r26sk_1 


(141) 


to  match  the  model.  Take  the  z-transforms,  setting  the  initial  condi¬ 
tion  vector  equal  to  zero: 


[l2  ”  ♦mi^m  "  rm^m 


(142) 


[iz  -  ♦aJXg  -  [fj  +  r2z~l]5s 


Given  a  control  law  of  the  form 


(143) 


"  cF«m  +  CBXS 


(144) 


Equation  143  becomes 


[Iz  -  *8  -  (Tl  +  r2z"1)cB]xg  -  [rL  +  r2z“  ]cF6m 


(145) 


The  relationships 


[rx  +  r2z-i]cF  -  rm 


(146) 


[rl  +  r2z”1]cB 


(147) 


cannot  be  identically  satisfied,  since  [Tj  +  does  not,  in 

general,  exists  Taking  a  mean-square  approach,  form  the  vector  E 
according  to  Equation  148, 


E 


trl 


+  ^z"1] 


(148) 


to  set  up  the  integrand  of  the  performand  index 


(149) 


The  gradient  of  E'(z)E(z)  with  respect  to  [CF(l/z)  |  CB(l/z)J  gives  the 
W-H  equation: 


-L 

z 


[r{  +  r'z] 


■v 

/ 

r 

F 

m 

CB 

w  -J 

1 

-e- 

3 

1 

-©- 

01 

I _ 

i 

* 


(150) 


This  verifies  the  W-H  equation  given  for  the  slewer  in  the  text. 
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APPENDIX  G 


-1 


woof,  c;‘cB  -  (rjr>)'1rj(*>  -  *a) 


The  relationship 


cFlcB  "  (rlrm)  irl(^m  ~  $s) 


%  Wr'i 


can  be  proven  by  using  the  properties  of  spectral  factorization. 


[V[  +  T2z){T1  +  T2z~l) 


have  the  factorization 


♦  ■ 


and  write 


^1*^1 [CF  I  cb1  -  ^  1  ,  ~  t rm  I  V  -  ♦s] 


The  symbolic  solution  is 

\-i  [rj  +  r2z] 


[CF  I  CB]  -  (fri1 


♦l* 


[^m  I  $m  “  ^s  ] 


But  and  Tj  +  V2z  have  no  poles  interior  to  the  unit  circle, 
fore,  only  the  pole  at  the  origin  can  contribute  to  the  partial 
expansion.  Thus, 


[cF  |  CB]  - 


-1 


4>li<0)  r£rm  |  ♦li(0)  rj(*m  -  ♦g) 


(151) 

Let 

(152) 

(153) 

(154) 

(155) 

There- 

fraction 

(156) 
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Since 


Cf1  -  (r{  r n,)"1  ♦u(O)*! 

then 

Cpts  "  (^lrm)-1  ♦l*(0)  4>1 4>1 1  <t»iiri ( 4»n 

or 

Cf^b  -  (rjr.)"1  rl(^  -  *8) 


APPENDIX  H 


PHYSICAL  REALIZATION  OF  THE  SLEVER  COUPLER 


This  appendix  reviews  the  conversion  of  a  sampled  sequence  to  a  con¬ 
tinuous  variable  using  a  zero-order  hold.  One  particular  realization  of 
the  ZOH  can  be  considered  to  be  an  Incremental  hold. 

A.  CONVENTIONAL  ZERO-ORDER  HOLD 

Consider  the  Impulse  sampling  of  a  time  signal  x  =  t,  with  a  trans¬ 
form  of 


x  *  1/s2 

Reconstruction,  via  a  ZOH  (Fig.  41a)  gives 


(157) 


Y 


1  -  e-sT 
s 


X(s) 


1  -  e~sT  Tz 
s  (z  -  1)2 


.  (158) 


Figure  41.  Zero-Order  Hold  Reconstruction 


119 


The  reconstruction,  using  Eq.  158,  Is  represented  schematically  In 
Fig.  41b. 


I  Figure  41b  represents  a  hardware  Implementation  using  a  clamping 

! 

>  circuit.  That  Is,  a  constant  level  of  voltage  (or  current)  is  held  for 

j  T  seconds  and  then  (physically)  removed.  The  control  law  then  (essen- 

I  tially  Instantaneously)  supplies  a  new  pulse,  which  is  clamped  (held) 

*  »eT 

,  for  the  next  T  second  interval.  The  data  hold  is  (1  -  e”  A)/s. 

i 


i 


Figure  42.  An  Incremental  Implementation 


B.  A  PARTIAL  SOFTWARE  IMPLEMENTATION 

Suppose  the  decision  is  made  to  relegate  the  transcendental  element 
of  the  clamping  circuit  to  software  (Fig.  42).  The  describing  equation 
18 


Y 


J. 

8" 


•  •  • 


(159) 


The  implication  is  (Fig.  42b)  that  each  T  second  an  increment  to  the 
control  action  is  added  to  the  signal  already  present.  The  data  hold  is 
1/s. 


:  ij  P  * 


The  preceding  viewpoints  transfer  readily  to  the  slewer  data  hold 


The  hardware  implementation  (Fig.  43)  interprets  the  output 


Figure  43*  A  Hardware  Slewer 

as  being  built  from  shifted  triangles  while  a  more  sof tware-oriented 
Implementation  might  retain  an  Integrator  and  a  ZOH  (In  hardware)  and 
relegate  the  remainder  to  software  (Fig.  44).  This  tends  (Fig.  44)  to 
have  the  appearance  of  an  Incremental  controller. 


Figure  44*  Partially  Software  Implemented  Slewer 


The  point  we  wish  to  make  at  this  juncture  Is  simple*  The  analyti¬ 
cal  model  used  In  the  synthesis  process  can  differ  from  the  physical 
Implementation;  however,  the  derived  control  law  properties  remain 
intact* 


APPENDIX  I 


ADDITIONAL  SLEWER  DATA  AND  TABULAR  LISTING 
OF  STEP  RESPONSES 


The  matrices  Involved  in  the  slewer  solution  are  $s»  rm,  r  ^ , 
and  I^.  These  are,  with  the  exception  of  Tj  and  Tj,  listed  in  Appen¬ 
dix  A*  Define 


R101 

r102 

r103 

R 1 2 1 

r122 

r123 

R104 

r105 

r106 

R124 

r125 

r126 

rl  - 

r107 

r108 

r109 

,  r2  - 

r127 

r128 

r129 

R1 10 

Rlll 

r112_ 

r130 

R 13 1 

r132  _ 

Data  for  Equation  162  are  listed  in  Table  7. 

TABLE  7 


DATA  FOR 

ri>  r2 

R101 

m 

-0.283742140 

R121 

Ml 

-0.197317335 

R102 

m 

0.115668810 

R122 

m 

0.218682553 

R103 

m 

0.102104000 

R123 

m 

0.095739417 

R104 

m 

0.008363868 

R124 

m 

0.014969735 

R105 

m 

-0.006567369 

R125 

m 

-0.013856154 

R106 

m 

-0.005131947 

R126 

m 

-0.008120770 

R107 

m 

0.332982000 

R127 

m 

0.277639404 

R108 

m 

0.294238000 

R128 

m 

0.300734839 

R109 

m 

2.287250000 

R129 

m 

2. 132551283 

R110 

m 

0.010283700 

R130 

m 

0.022925586 

Rlll 

m 

0.009357000 

R131 

m 

0.020065493 

R112 

m 

0.076482000 

R132 

m 

0.152324084 

TABLE  8.  SLEWER  STEP  RESPONSES 


Ml'  0.00080W30  0.845737063 
M2'  I.0MWMM  1.0881 176V 
M3'  O.000000000  -8.813466723 
M4»  0.000000000  -4.013747898 


Ml'  0.867496861  0.040147823 
802*  0.801243138  8.088368214 
883'  8.023416892  8.003729333 
804»  -0.088983360  -8.008884114 


0.860434968 

0.888621334 

-0.802199976 

-8.000728331 


0.060677113  0.868783469 
0.088322119  0.803646123 
-0.000103933  0.000021235 
-8.081329384  -0.801263243 


RI3'  1.333436216  4.696242196 
114'  1.667517894  0.665216749 
015'  -0.458339559  -0.172334136 
0.  3. 


>13=  1.147232834 

8.394139300 

1.837649412 

8.972434251 

8.977715595 

>14=  1.229792917 

0.S94S26302 

1.107724311 

0.337333531 

1.834175433 

>15=  -0.333396952 

-0.245963690 

-0.297762217  -0.271497306 

-0.277453132 

18. 

15. 

20. 

25. 

38. 

001'  4.817379333  0.870832452 
M2'  8.884292115  0.806428972 
>03'  -0.030708233  8.088608122 
>84'  -0.083952695  -0.016446332 


>81= 

0.865105370 

0.068130317 

8.860366755 

8.868364974 

0.066629634 

>82= 

8.000365533, 

0.060254143 

0.088348312 

0.806656734 

0.800763639 

>93= 

0.823563335 

0.802019039 

-0.001553394 

-0.600215331 

0.880233323 

>84' 

-8.006496713 

-0.808429312 

-0.880928712 

-8.081313972 

-0.881296336 

>13'  0.562121392  1.160622941 
>14'  0.343256696  1.249593669 
>15*  -0.139155544  -0.334343938' 
I.  6. 


>13=  0.34*837595  1.073152247 
>14=  0.335675474  1.155334863 
>15=  -0.229132429  -0.311752351 

11.  14* 


0.944918031  8.996671833  1.801599991 

0.954267118  1 . 057959539  1 . 023927932 

-0.262222748  -0.233994441  -0.231518494 

21.  26.  31. 


Ml'  0.034367119  0.871133194 
>02=  0.007522752  0.085149356 
>83=  -0.039460099  0.013653962' 
M4*  -0.006202290  -0.015661312 


>01=  0.063568400  0.059793991 
>02=  0.990249139,  0.800534926 
>03=  0.018731427'  0.000344499 
>84=  -0.004339853  -0.888316389 


0.960713305 

8.809672744 

-0.081594194 

-0.081964332 


9.860676733 

0.060795131 

0.099237276 

-0.891345226 


8.068761370 

8.009652158 

-8.000934765 

-0.091241732 


>13=  1.243711064  0.754954354 
>14=  t. 351203724  0.723347837, 
>15=  -1.360440473  -0. 194884193' 
2.  7. 


>13*  1.126732450  >.*12353795  1.021350003  8.933746710 

>14=  1.299111486  8.316148319  1.93S715737  1.901233371 

>15*  -0.327336399  -8.251714261  -0.232333530  -0.275369146 
12.  17.  22.  27. 


8. 978210238 
1.924476647 
0.274357915 
32. 


Ml*  9.043953173  0.071952527 
M2*  0.009824367  0.903491572 
M3*  -0.033337352  0.921176209 
M4*  -9.910243044  -0.013946238 


Ml*  0.061817274  8.868121379  0.869592219  0-060329573  0.060624713 
>02=  0.090327604  0.808479761  0.000847060  0.000848621  8.008764762 
003*  0.814521491  -0.001439832  -8.080731215  0.008006470  0.900145433 
>04=  -9.092722226  -0.080369937  -0.001287346  -0.001293267  -0.901275251 


>13*  1.641603235  1.160535731 
114'  0.613504065  1.241354598 
>13'  -0.155571322  -0.335335381 
3.  8. 


>13=  9.872908715  1.856326091  0.959443596  0.986503168 

>14=  0.869407833  1.129906334  8,971533820  1.045306317 

>15*  -0.233930097  -0.304193961  -8.267835943  -0.230461021 
13.  18.,  23.  23. 


1.003447031 

1.032901599 

-0.283397802 

33. 


Ml'  0.859004422  0.869200935 
M2'  0.003774122  0.082435093 
M3*  -0.023722649  9.935330934 
104'  -0.813674109  -0.011538742 


>81=  0.061073213  0.060028565  0.060347403  8.068649301 

>82=  0.008072326  0.000761434  8.800672025  0.000777393 

>03*  0.899217262  -0.081629300  .9.999793155  0.000301236 

>04=  -8.001549529  -0.800542489  -0.081257643  -0.00131325* 


0.068743993 

0.000660339 

-0.000032491 

-8.081232690 


>13=  1.184939332  0.886323102 
>14*  1.271957936  0.73’!55t2l 
>U»  -0.333276641  -0.214575351 
4.  9. 


>13*  1.102411303  0.929134323  I.083353222  0.993421705 

>14=  1 . 1 S273525J  8.935729663  1.97221  1418  1.013446629 

>15*  -0.320924325  -0.257063637  -0.287333363  -0.273695163 
14.  1*.  24.  2». 


0.963917495 

1.8161484*7 

-0.272657047 

34. 


8.868630303 

0.000762008 

8.088835864 

-0.001272412 

1.814113340 

1.048527945 

-8.235334349 

35. 


0.060744902 

8.000665051 

-0.880096414 

-0.8012331:6 

0.953749276 

1.009O33263 

-0.270829096 

36. 


0.060641853 

0.080753592 

0.800057272 

-0.001275385 

1.013726043 

1.046942329 

-0.287523465 

37. 


0.060743823 

0.000674231 

-0.080039615 

-0.081236393 

0.954604253 

1.003204223 

-0.269348393 

33. 


8.060653519 

0.000754735 

0.000045019 

-0.001275361 

1. 0223:9042 
1.052253213 
-8.233S49013 
39. 
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